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Classification of symmetric periodic trajectories in ellipsoidal billiards 
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We classify nonsingular symmetric periodic trajectories (SPTs) of billiards inside ellipsoids of IR"+^ without 
any symmetry of revolution. SPTs are defined as periodic trajectories passing through some symmetry set. 
We prove that there are exactly 2^"'(2"+^ — 1) classes of such trajectories. We have implemented an algorithm 
to find minimal SPTs of each of the 12 classes in the 2D case {M?) and each of the 112 classes in the 3D case 
(K^). They have periods 3, 4 or 6 in the 2D case; and 4, 5, 6, 8 or 10 in the 3D case. We display a selection 
of 3D minimal SPTs. Some of them have properties that cannot take place in the 2D case. 



Smooth convex billiards are a paradigm of con- 
servative dynamics, in which a particle collides 
with a fixed closed smooth convex hypersurface 
of M""*" . They provide examples of different dy- 
namics: integrable, mostly regular, chaotic, etc. 
In this paper we tackle the integrable situation. 
Concretely, we find and classify symmetric peri- 
odic trajectories (SPTs) inside ellipsoids of R""^ . 
PTs show different dynamics, which describe how 
they fold in M""*""*^. STs present symmetry with re- 
spect to a coordinate subspace of M""*" . Dynamics 
and symmetry are precisely the main aspects we 
consider in our classification of SPTs. We estab- 
lish 112 classes of SPTs in the 3D case, and we 
find a representative of each class with the small- 
est possible period. Those minimal SPTs have 
periods 4, 5, 6, 8, or 10. We depict a selection 
of minimal 3D SPTs. Some of them have proper- 
ties that cannot take place in the 2D case. SPTs 
are preserved under symmetric deformations of 
the ellipsoid. In a future paper we plan to study 
their bifurcations and the transition between sta- 
bility and instability under such deformations. 



I. INTRODUCTION 

Smooth convex billiards, in which a point particle 
moves uniformly until it undergoes abrupt elastic col- 
lisions with a smooth strictly convex hypersurface Q C 
M"+^, are one of the typical examples of conservative sys- 
tems with both regular and chaotic dynamics. The 2D 
case (n = 1) was introduced by Birkhoff^ and there is 
an extensive literature about it. The high-dimensional 
case remains much less studied, although we highlight 
two known results. First, there exists a Nekhoroshev-like 
theorem for billiards.^ Second, there are some general 
lower bounds on the number of periodic billiard trajecto- 
ries. '^^^ Of course, these lower bounds are not informative 
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for integrable systems, where periodic trajectories are or- 
ganized in continuous families. 

Ellipsoids are the only known examples of smooth 
strictly convex closed hypersurfaces of the Euclidean 
space M"'+^ that give rise to integrable billiards. There is 
a larger class of integrable billiards, called Liouville bil- 
liard tables, but in none of the new integrable billiards 
the point particle follows straight lines between two con- 
secutive impacts; instead it moves on the geodesies of 
some smooth compact connected Riemannian manifold 
with boundary. ^^^ 

In fact, an old conjecture atributed to Birkhoff states: 
ellipsoids are the only smooth strictly convex closed hy- 
persurfaces of M""*"^ that give rise to integrable billiards. 
The conjecture remains still unproven, but Berger^" and 
Gruber^^ proved that only ellipsoids have smooth caus- 
tics when n > 2, which is a powerful argument for the 
conjecture. A smooth caustic is a smooth hypersurface 
with the property that a billiard trajectory, once tangent 
to it, stays tangent after every reflection. 

Periodic trajectories (PTs) are the most distinctive tra- 
jectories, so their study is the first task. There exist two 
remarkable results about PTs inside ellipsoids: Poncelet 
theorem and Cayley conditions. Poncelet"'^^ showed that 
if a billiard trajectory inside an ellipse is periodic, then all 
the trajectories sharing its conic caustic are also periodic. 
Cayley^'^ gave some algebraic conditions for determining 
conic caustics whose trajectories are periodic. Later on, 
both results were generalized to any dimension. ^"^'^^ 

Almost all papers on billiards inside ellipsoids follow 
an algebro-geometric approach to the problem. ^^^^'^ A 
more dynamical approach to billiards inside ellipsoids has 
recently been considered in Refs. 24-27. We follow that 
approach, in which the description of the dynamics is 
more important than the obtaining of explicit formulae. 
Indeed, this paper is a natural continuation of the study 
on (minimal) periodic trajectories carried out in Ref. 27. 
In the framework of smooth convex billiards the minimal 
period is always two. Nevertheless, since all two-periodic 
billiard trajectories inside ellipsoids are singular — in the 
sense that some of their caustics are singular elements 
of a family of confocal quadrics — , two questions arise: 
which is the minimal period among nonsingular billiard 
trajectories? which ellipsoids display such trajectories? 
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FIG. 1. Two SPTs of period 6, but with only 5 and 4 different impact points. Lines in green and yellow represent the 
intersections of the ellipsoid with Hl-caustics and H2-caustics, respectively. The E-caustic is also shown in the EH2 figure. 



The authors settled both questions in Ref. 27. 

We say that an ellipsoid of M"+^ whose n + 1 axis 
have different lengths (i.e., an ellipsoid without any sym- 
metry of revolution) is nondegenerate. Our current goal 
is to classify nonsingular symmetric periodic trajectories 
(SPTs) inside nondegenerate ellipsoids. We also look for 
minimal SPTs — representatives of each class with the 
smallest possible period. 

The relevance of the SPTs relies in three facts. First, 
they persist under small enough symmetric perturbations 
of a nondegenerate ellipsoid. Second, each perturbed 
SPT can be computed by solving a nonlinear equation for 
which the corresponding unperturbed SPT is a good ini- 
tial approximation. Third, this nonlinear equation is n- 
dimensional, not 2n-dimensional. That is, the dimension 
of the problem can be halved. These properties greatly 
simplify the numerical computation of SPTs inside sym- 
metrically perturbed ellipsoids, which is our next aim. 
To be more precise, we plan to study the stability and 
bifurcations of perturbed SPTs, as a way to gain some 
insight into the transition from integrability to chaos in 
high-dimensional smooth convex billiard problems. 

Roughly speaking, a symmetric trajectory (ST) inside 



Q 



(a;i,...,x„+i) e 



pn+l 



ai 



a-n+i 



is a billiard trajectory invariant under a orthogonal re- 
flection, with respect to some of the 2"+^ — 1 coordinate 
proper subspaces of the Euclidean space M"+^. For in- 
stance, the left trajectory in Fig. 1 is invariant under the 
xi-axial and xiXs-specular reflections, whereas the right 
one is invariant under the X2-axial reflection. There are 
twice as many "billiard symmetries" as coordinate proper 
subspaces of M""*"^. This is easily understood by means 
of the following example. We can associate two different 
kinds of SPTs to any coordinate axis of M""*"^; namely, 
the ones with an impact point on the axis and the ones 
with a segment parallel to the axis. Note that the left 
SPT in Fig. 1 has an impact point on the xi-axis — and 



so, it is invariant under the xi-axial reflection — , but also 
has a segment parallel to the X2-axis — and so, it is in- 
variant under the Xi^s-specular reflection. We will check 
that this kind of double symmetry is the norm among 
SPTs, not the exception. 

In view of the above arguments it is easy to guess that 
there is a rich casuistry of SPTs inside nondegenerate 
ellipsoids of IR"+^, so much so that there is no simple way 
to describe them. For simplicity, let us briefly consider 
the 3D case: n — 2. Then, any billiard trajectory inside 
Q C M'^ has two caustics of the form 



Q\ = { {xi,X2,X3) £ 



3 

E 



Let us assume that < oi < a2 < as. We restrict our 
attention to nonsingular trajectories. That is, trajecto- 
ries whose caustics are ellipsoids: < A < ai; 1-sheet 
hyperboloids: oi < A < 02; or 2-sheet hyperboloids: 
02 < A < 03. The singular values A € {ai,a2,a3}, are 
discarded. It is known that there are only four types of 
couples of nonsingular caustics: EHl, HlHl, EH2, and 
H1H2. The notation is self-explanatory. It is also known 
that each caustic type gives rise to a different billiard dy- 
namics. For instance, only trajectories with EH2-caustics 
rotate around the xa-axis. The left trajectory in Fig. 1 
is a sample. 

Coming back to the general case, its billiard trajecto- 
ries present 2" different caustic types, which describe how 
trajectories fold in M"+^; see Ref. 27. Then, we classify 
SPTs inside nondegenerate ellipsoids of ]R"+^ by their 
caustic types and their symmetries. That is, two SPTs 
with the same symmetries, but different caustic type, are 
in different classes. On the one hand, not every symmetry 
takes place among SPTs of a given caustic type, but only 
2"+^. On the other hand, all SPTs are, in some sense, 
doubly symmetric, which rise to 2" (2"+^ — l) classes we 
may encounter with a given caustic type. Consequently, 
we know that there are at most 2^" (2"+-'^ — l) classes of 
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SPTs. Finally, by means of an argument involving wind- 
ing numbers, we establish the existence of SPTs for all 
those classes. 

That concludes the analytical part of the work. We 
have also implemented a numerical algorithm to find min- 
imal SPTs of each of the 12 classes in the 2D case, and 
each of the 112 classes in the 3D case. They have periods 
3, 4 or 6 in the 2D case; and 4, 5, 6, 8 or 10 in the 3D 
case. We depict all 2D minimal SPTs in Table V, and 
just a gallery (for brevity of exposition) of 3D minimal 
SPTs in Sec. VI. 

For instance, we show in Fig. 1 two minimal SPTs of 
period 6, but with only 5 and 4 different impact points. 
Only the second phenomena can take place in the 2D 
case. This has to do with the fact that, from a generic 
point on the ellipsoid, we can trace four lines tangent to 
a fixed couple of caustics. On the contrary, we can trace 
just two in the 2D case. 

Some of the ideas used in this paper were first in- 
troduced by Kook and Meiss,^® to classify symmet- 
ric periodic orbits of some standard-like reversible 2n- 
dimcnsional symplectic maps. The technical details of 
our problem are harder — there are more symmetries, 
caustic types play a role, and the classification does 
not depend only on the evenness or oddness of some 
integers — , but the main arguments do not change. 

Some of our billiard SPTs could be considered as dis- 
crete versions of the closed geodesies on triaxial ellipsoids 
found in Ref. 29. For instance, it is interesting to com- 
pare the SPT on the left side of Fig. 1 with the symmetric 
closed geodesic shown in Fig. 4 of that paper. Of course, 
the SPT on the right side of Fig. 1 have no continu- 
ous version, because there are no geodesies with return 
points. 

We complete this introduction with a note on the or- 
ganization of the article. We review the classical theory 
of reversible maps in Sec. II following Ref. 30. Next, 
we specialize that theory to billiards inside symmetric 
hypersurfaces. In Sec. IV we review briefly some well- 
known results about billiards inside ellipsoids, in order 
to fix notations that will be used along the rest of the 
paper. Billiards inside ellipses of M^ and inside triaxial 
ellipsoids of M*^, are thoroughly studied in Sec. V and 
Sec. VI, respectively. Billiards inside nondegenerate el- 
lipsoids of M""*"^ are revisited in Sec. VII. Perspectives 
and conclusions are drawn in Sec. VIII. 



II. REVERSIBLE MAPS AND SYMMETRIC ORBITS 

In this section we state the concept of symmetric pe- 
riodic orbit for a map. A map is reversible if each or- 
bit is related to its time reverse orbit by a symmetry 
transformation. Reversible maps can be characterized as 
maps that factorize as the composition of two involutions. 
Symmetric orbits must have points on certain symmetry 
sets; namely, the fixed sets of the involutions (reversors) 
which factorize the original map. A reversible map can 



have different factorizations and, therefore, its symmetric 
orbits can be classified according to the symmetry sets 
they intersect. 

Let / : M ^ M be a diffeomorphism on a manifold 

M. 

Definition 1 The map f is symmetric when there exists 
an involution s : M ^ M (i.e., s o s ~ Id), such that 
f o s — so/. Then, s is called a symmetry of the map f . 

Definition 2 The map f is reversible when there exists 
an involution r : M —> M such that for = rof^^. Then, 
r is called a reversor of the map f (i.e., f is r -reversible). 

Remark 1 The composition of a symmetry and a rever- 
sor (provided both commute) is another reversor. 

Definition 3 If r is a reversor, we denote by Fix(r) = 
{m Cz M : r{m) — m} its set of fixed points, which is 
called the symmetry set of the reversor. 

Definition 4 An orbit of the map f is a sequence O = 
{mj)j^z such that mj — /(?7ij_i) = /■'(ttio). 

Definition 5 An orbit O of a r-reversible map f is called 
r-symmetric when r{0) ~ O . 

The following characterization of reversible maps goes 
back to G. D. Birkhoff. 

Lemma 1 A map f is reversible if and only if it can be 
factorized as the composition of two involutions, in which 
case both of them are reversors of f. 

Proof Let us assume that / is f-reversible. Then f = 
f o f ~ r o f^^ is another reversor, because: 



f O f = f I 



r^ = fof- 



and 



mf~fof^rof o f o f = r = Id. 

Therefore, the map / = / o f^ = f o f is the composition 
of two involutions. 

On the other hand, if / = for and f'^=r'^= Id, then: 

• /ofo/ = fof^ofof = f^of = f, and 



• /ofo/ = fofofof = fof=f, 

so both involutions r and f are reversors of /. 



D 



The factorization of a reversible symmetric map as a 
composition of two involutions is not unique. From any 
given factorization / = f o f , we can construct infinitely 
many more, namely, f ^ Vk o fk, where fk — f^'^^ o f — 
f^ of = fof~^ and fi^ — f^of for any k E Z. Neverthe- 
less, these new factorizations do not provide new sym- 
metric orbits, because f{0) ~ r{0) = fk{0) = rk{0), 
i.e., an orbit is invariant under f (or f) if and only if it 
is invariant under all the reversors fj, and fj,. On the 
contrary, the existence of a symmetry s commuting with 
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the reversors f and f is more promising, because then 
/ = (s o f) o (s o f) is a new factorization, which could 
give rise to new kinds of symmetric orbits. The more such 
symmetries exist, the more kinds of symmetric orbits we 
can try to find. 

Given a factorization / = f o f , we say that f and f 
are associated reversors, and Fix(f ) and Fix(f) are asso- 
ciated symmetry sets. The importance of these concepts 
is clarified in the following characterization of symmetric 
orbits, which can be found in Ref. 30. 

Theorem 1 Let f = for be any factorization of a rever- 
sible map as the composing of two involutions. Then an 
orbit of this map is: 

a) f -symmetric if and only if it is f -symmetric. 

b) f -symmetric if and only if it has at least one point on 
Fix(f)UFix(f), in which case it has no more than two 
points on Fix(f) U Fix(f). 

c) f -symmetric and periodic if and only if it has exactly 
two points on Fix(f) U Fix(f), in which case it has a 
point on each .symmetry set if and only if it has odd 
period. In particular, 

i. An orbit is f -symmetric with period 2k if and only 
if it has a point m^ such that either m^ G Fix(f) n 
/ Fix(f) or TO* G Fix(f) n / Fix(f), in which 
case f*'{m^,) is on the same symmetry set as to.*. 

ii. An orbit is f -symmetric orbit with period 2k + 1 
if and only if it has a point m^ such that to,* e 
Fix(f) n /'^Fix(f), in which case f'''^^{m^) G 
Fix(f). 

Proof This is an old result, so we skip most of the 
details. We focus in those aspects that are more relevant 
for our goals. 

a) It is immediate as f (O) = f (O). 

b) To begin with, let us assume that an orbit O ~ 
{mj)jQZ is f-symmetric, so there exists some / G Z 
such that mi — f(?Tio). Then mi^i — f{mi) = f{mQ), 
because f — f o f. We distinguish two cases: 

• If Z = 2/c is even then ruk G Fix(f) as: 

mk = .r''{m2k) = .r''(r(mn)) = f{mk). 

9 li I = 2k + 1 is odd then nik+i G Fix(f) as: 

TOfc+i = /-('=+in"^2fe+2) - f-^'+'Hri^o)) = H^k+l)- 

On the other hand, let us assume that O has a point, 
say mk, in one of those fixed sets. Then there exist 
two possibilities. In the case to^ G Fix(f), Vj G Z: 

rrik+j = f^{mk) = f'{f{mk)) = f{f^'{mk)) = f{mk-j). 



In the case ruk G Fix(f), from f = f ^ o f, \fj G Z: 

mk+j-i = .P^^imk) = f'^^{f{mk)) 

= /^(f(mfe)) = f{f^^{mk)) = f{mk-j). 

Hence, O turns out to be f-symmetric in both cases. 

c) Next, we study the symmetric periodic orbits. First, 
if O has two points, say TOq and m^, on Fix(f), then 
O is 2A:-periodic, because: 

m2k = .f{mk) ^ I^{f{mk)) 

^f{f^^{mk)) == f(TOo) = TOq. 

Second, if O has one point on each symmetry set, say 
Too G Fix(f) and m^^i G Fix(f), then O is {2k + 1)- 
periodic, because: 

TO2fe+l = /''(TOfe+l) = /*=(f(TOfe + l)) 

= /'=+l(f(TOfe+l))=f(r('=+l)(TOfc+l)) 

— f{mf)) = JTlg. n 

Therefore, the computation of symmetric orbits with 
period 2k can be reduced to the computation of points 
on symmetry sets, which are mapped onto the same sym- 
metry set after k iterations of the map. On the other 
hand, the computation of symmetric orbits with period 
2fc + 1 can be reduced to the computation of points on 
symmetry sets, which are mapped onto their associated 
symmetry sets after k -\- 1 iterations of the map. 



III. BILLIARDS INSIDE SYMMETRIC 
HYPERSURFACES OF E"+^ 

Let Q be a (strictly) convex smooth hypersurface of 
jjn+i rpj^g billiard motion inside Q can be modelled by 
means of a diffeomorphism / defined on the phase space 

M = {{q,p) G Q X §" : p is directed outward Q a.t q} . 

We define the billiard map f : M —>■ M, f{q,p) — (q' ,p'), 
as follows. The new velocity p' is the reflection of p, with 
respect to the tangent hyperplanc TqQ. That is, 

p' =Pt -Pn ==p-2p„, 

where pt and Pn are the tangent and normal components 
at q of the old velocity: 

P == Pt+Pn, Pt e TqQ, Pn G NqQ . (1) 

Let q+ {p') be the line through q with direction p' . Then 
the new impact point q' is determined by the condition 

Qn{q+{py) = {q,q'}. 

That is, q' is the intersection of the ray {q + fip' : /i > 0} 
with the hypersurface Q. This intersection is unique and 
transverse by convexity. 
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Definition 6 A billiard orbit is a sequence of points 
{nij} -^^ C M such that rrij — /(toj„i) — /-'(ttio). // 
the points nij ~ {qj,pj) form a billiard orbit, then the 
sequence of impact points {qj} 



jez. 



C 



is a billiard (or 
impact) configuration, whose joining by polygonal lines 
form a billiard trajectory; the sequence of outward veloc- 
ities {pj} ^2 C §" is a velocity configuration. 

The distinction between orbits and trajectories is clear. 
We refer to orbits when we are working in the phase space 
M, whereas we refer to trajectories when we are draw- 



ing m 



pn+l 



There exists an one-to-one correspondence 
between billiard orbits and billiard configurations. The 
velocities of the billiard orbit that correspond to the bil- 
liard configuration {qj} ^^ ^'"^• 



Pj 



<1] - gj-i 



(2) 



From now on we shall assume that the hypersurface Q 
is symmetric with respect to all coordinate hyperplanes 
of IR"+-'^. Then the billiard map admits 2"+^ factoriza- 
tions as a composition of two involutions. We need the 
following notations in order to describe them. 

Let S be the set made up of the reflections — that is, 
involutive linear transformations — with respect to the 
2"+^ coordinate subspaces of IR"+^. We represent its 
elements as diagonal matrices whose diagonal entries are 
equal to 1 or —1, so 

S = {(7 = diag(CTi, . . . ,cr„+i) : cfj e {-1, IjforaUj} . 

Given {q,p) G M, there exists a unique point q such that 

Qniq+(p)) = {q,q}. 

Thus, q denotes the previous impact point. Finally, let 
jj be the reflection of the velocity p with respect to the 
normal line NqQ. That is, 

P = -p' = Pn - Pt = P - 2pt, 

where pt and Pn were defined in the decomposition (1). 

In what follows, symbols g', p', q, and p have the mean- 
ing given in the previous paragraphs. We emphasize that 
they make sense only after both impact point q, and uni- 
tary outer velocity p, are given. Next we describe the 
2"+^ factorizations of the billiard map. 

Proposition 1 Let f : M ^ M, f {q,p) = {q',p'), be 
the billiard map inside a closed convex symmetric hyper- 
surface Q C M"^"'^. Let f,f:M^M be the maps 

f{q,p)^ {q,p), f{q,p) = {q,-p). 
Let Sa,f„,f„ : M ^f M be the maps 

Sa{q,P) = {<yq,crp), fa^S^of, f^^S^of 

for any reflection cr G S. Then: 



a) f ~ f o f, and f 2 _ ^2 _ j^-^ 

b) f o Sa ^ s^o f, si ^ Id, fos^^s^of, and f o s^, ^ 
Sa o r. 

c) f ~ f^ o fcr, and r1 = f1 = Id. 
Proof 

a) r{f{q,p)) = f{q,p) = f{q,-p') = {q' ,p'). Besides, 
P{q,p) = f{q,-p) = {q,p), as Q n {q + {-p)) = 
{q,q}. Finally, f'^iq,p) =r{q,p) = {q,p). 

b) Clearly s^ is an involution. Since Q is symmetric with 
respect to the coordinate hyperplanes, it turns out 
that if Pt and p^ are the tangent and normal compo- 
nents of a velocity p at an impact point q, then apt 
and apn are the tangent and normal components of 
ap at the impact point aq. On the other hand, if we 
write q' = q + fi{q,p') p' , then n{crq,ap') = fJ.{q,p'), 
again by the symmetry of Q. Hence, f o s„ = s^ o f. 
The proof oi f o Sa- = s^. of and f o s^ — s„ of follows 
the same lines. 



c) It is a direct consequence of Remark 1. 
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No symmetry has been required to obtain the factor- 
ization / = for. Therefore, all convex hypersurfaces 
Q C M"+^ give rise to reversible billiard maps, although 
symmetric hypersurfaces have much more factorizations. 

We introduce the acronyms SO, ST, SPO, and SPT 
for symmetric orbit, symmetric trajectory, symmetric pe- 
riodic orbit, and symmetric periodic trajectory, respec- 
tively. If r is any reversor of the billiard map, we can 
deal with r-SOs, r-STs, r-SPOs, and r-SPTs. 

Once these 2"+^ factorizations f = fcr ° fa have been 
found, we describe their symmetry sets. In the next 
proposition we prove that only two symmetry sets arc 
empty. We also provide an explicit geometric description 
of the 2"+^ — 2 nonempty symmetry sets where one can 
look for SOs. As a consequence, we shall see that these 
symmetry sets are mutually exclusive — that is, they do 
not intersect — but at some very specific points which 
are described in detail. As before, some notations are 
required. 

We recall that given a reflection defined on an Eu- 
clidean space, we can decompose the space as the or- 
thogonal sum of the eigenspaces of eigenvalues — 1 and 1 . 
This is called the spectral decomposition of the reflection. 
Let M"+^ = V^ _L V~ be the spectral decomposition 
associated to any reflection cr G S. That is. 



V^ = {p G 



pri+l 



ap = ±p}. 



For instance, if (xi, X2, x^) are the Cartesian coordinates 
in R^ and a = diag (—1, 1, 1), then V^ is the a;2a;3-plane, 
and V-T is the xi-axis. We also introduce the sections 



Q, = Qnv+ 
p„ = §"nTC 



= {qeQ ■.aq = q}, 
= {p G §" : crp = -p}. 
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Given two linear varieties Vi,F2 C K""'""'^, the symbol 
Vi _L* V2 means that they have a orthogonal intersection. 
Finally, we recall that a line in M"+^ is a chord of Q when 
it intersects orthogonally Q at two different points. 

Proposition 2 The symmetry sets of the reversors f^ 
and r„, cr G S, are: 

Fix(f^) = {{q,p) eM -.qeQa andp G NgQ„} , 
Fix(f^) = {{q,p) e M ■.<7qeq+ (p) andp G P^} 

^{{q,p)eM:{q+(p))UV+}. 

Only the reversors r-u and fjd have empty symmetry 
sets. Moreover, if a point {q,p) G M belongs simultane- 
ously to two different symmetry sets, then the line q+{p): 

• Is contained in some coordinate hyperplane of M."'^^ , 
in which case so is its whole billiard trajectory; or 

• Is a chord of the hypersurface Q through the origin, in 
which case its billiard trajectory is 2-periodic. 

Proof To begin with, we deduce from the definitions of 
reversors f„ and f^, and symmetries s„ that 



(gf,p) G Fix(f^) 
{q,p) G Fix(f„) 



Sa{q,v) =r(q,p) 

<^q = q] ^^^^ fq e 

ap = p\ 1 crp = 

Sa{q,p) =r{q,p) 



aq = q 
up = — p 



crq e q + 



(P) 



To understand the condition ap — p, we compare the 
spectral decomposition M"+^ = Vj' _L V^T with the spec- 
tral decomposition R"+^ ~ TgR"+'^ = 
ated to the reflection p 1-^ p. We note that 



TgW^+' = V+ _L V- associ- 



V- - {p G r,M"+i ■.p=-p} 
V+ - {p G r,M"+i : p = p} 



If 9 G Qcr, then 
NgQ, = NgQ ± V- = % 



.9 



V- 



TqQa = TqQnv+ = v-nv+, 



N„Q. 



V- c V-. 



The equality on N^Q^ is due to {AOB)^ = A^LB^. 
The inclusions follow from V^ being a symmetry sub- 
space of Q. Those relations imply that 



rgM"+i 



N„ 



- {v+ n v+) ± {V- nv-)i. {V- n v+). 

Therefore, if we fix any q G (5<j, then crp ~ p ^ p Cz 

{v+ n v+) ± (v- n V-) = v+ ± v- = n,q„, so 

Fix(f,) = {{q,p) eM:qeQ„,pe NgQ,} . 

Next, we study the other symmetry sets for a 7^ Id. If 
p (z Pa, then p _L Vj^ and p || V^ , which implies that 
aq e q + (p) 4^ q + (p) ±* V+. Therefore, 

Fix(f,) ^{{q,p)eM:q+ (p) U V+} . 



Clearly, Q^ is empty if and only if cr = —Id, and P^ is 
empty if and only if ct = Id. 

Finally, let us check that if the point {q,p) G M belongs 
to two different symmetry sets, then q + (p) is a chord of 
Q or is contained in some coordinate hyperplane like 



Hj = {q= {xi,.. .,x„+i) G 



= 0} . (3) 



First, if (g,p) G Fix(fcr) n Fix(fT-), then using the rela- 
tions written at the beginning of this proof, we get 



aq 



Tq, 



up ~ p = Tp. 



If (T 7^ r, then (jj 7^ Tj for some j, so that q ^ Hj, p\\ Hj, 
and q+ {p) d Hj . 

Second, if {q,p) G Fix(fo-) nFix(fT), then using again 
the same relations, we get that 



aq = q, 



ap = p = — p 



Tq = q, 



rp = —p. 



If (T 7^ Id, then aj = —1 for some j, so that q,q (z Hj, 
p\\Hj, and q + (p) C Hj. If cr = Id but r 7^ —Id, then 
Tj — 1 for some j, so that p' — —p G NqQ, p\\ Hj, and 
p' \\ Hj. Thus, pn 1 1 Hj and NqQ C Hj , which implies that 
q G Hj, because Q is symmetric, smooth, and (strictly) 
convex. Therefore, q + (p) C Hj. If cr = Id and t = —Id, 
then p' = — p G NgQ and q = —q, so that the line q + (p) 
intersects orthogonally Q at q and —q. 
Third, if {q,p) G Fix {f^) H Fix (f^), then 



aq 



Tq, 



ap 



-p = TP, 



and we apply the same reasoning as in the first case. D 

Hence, only very specific billiard orbits — the ones con- 
tained in a coordinate hyperplane or 2-periodic — can 
have a point in the intersection of two symmetry sets. 
For instance, it turns out that any 2-periodic billiard or- 
bit inside a nondegenerate ellipsoid is contained in more 
than half of the 2"+^ — 2 nonempty symmetry sets. 

Nevertheless, there exists other billiard orbits that in- 
tersect different symmetry sets at different points. In- 
deed, all SPOs with odd period have points on two asso- 
ciated symmetry sets, whereas all SPOs with even period 
have two points on some symmetry set; see item c) of 
Theorem 1. However, a especial situation arises when Q 
is an ellipsoid: most of its SPOs with even period have 
exactly four — instead of the expected two — points in 
the symmetry sets described in Proposition 2, in which 
case they must intersect two different symmetry sets (see 
item b) of Theorem 1). 

That motivates the following definition. 

Definition 7 An SPO inside a symmetric hypersurface 
Q C M"'*'"'^, is a doubly SPO when it intersects two dif- 
ferent symmetry sets. 

Remark 2 Let O be a f^-SO or, equivalently, a f^-SO 
(see Theorem 1). By definition, O, as a subset of the 
phase space M , is invariant under the reversors 

fa{q,p) ^ {crq,crp), fa{q,p) == {aq,-ap). 
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In particular, the billiard configuration associated to 
O, viewed as a subset of the configuration space Q, is in- 
variant under the map a\Q : g i— )■ aq, whereas its velocity 
configuration, viewed as a subset of the velocity space §", 
is invariant under the map — i7|sn : p i~>- —ap. This mo- 
tivates the following definitions for the 2D and 3D cases. 
A billiard configuration inside a symmetric curve/surface 
is called central, axial or specular when it is symmetric 
with respect to the origin, some axis of coordinates or 
some plane of coordinates, respectively. Similar defini- 
tions apply to velocity configurations. 



IV. BILLIARDS INSIDE ELLIPSOIDS OF 



pn + l 



The billiard dynamics inside ellipsoids has several im- 
portant properties. For instance, it is completely inte- 
grable in the sense of Liouville. We present some results 
about such billiards. First, we list the classical ones, 
which can be found in the monographs Refs. 31-33. Next, 
we describe a dual property discovered in Rcfs. 14 and 16. 
Finally, we detail the behaviour of elliptic coordinates of 
billiard trajectories inside ellipsoids given in Rcf. 22. 



A. Caustics and elliptic coordinates 

The following results go back to Jacobi, Chasles, Pon- 
celet, and Darboux. We consider a billiard inside the 
ellipsoid 



Q={q^ 



pn+l 



{Dq,q)^l], 



(4) 



where D^^ =^ diag (ai, . . . ,a„+i) is a diagonal matrix 
such that < ai < • • • < a„+i. The degenerate cases in 
which the ellipsoid has some symmetry of revolution arc 
not considered here. This ellipsoid is an clement of the 
family of confocal quadrics given by 



<9m = {^e 



yri + l 



{D^q,q) = l], ne 



where Df^ = \^^^ ^ mW) . The meaning of Q^ in the 
singular cases fx € {ai, . . . ,a„+i} must be clarified. If 
fi = Uj, we define it as the n-dimensional hyperplane (3). 

Theorem 2 Let Q be the nondegenerate ellipsoid (4). 

a) Any generic point q G M""*"^ belongs to exactly n + 1 
distinct nonsingular quadrics Qf^^ , ■ ■ ■ , Qp.^ such that 

^0 < fll < Ml < 0^2 < • • • < On < Mn < O-n+l- BcsidcS, 

those n + 1 quadrics are mutually orthogonal at q. 

b) Any generic line £ C IR"'*'"'^ is tangent to exactly n dis- 
tinct nonsingular confocal quadrics Q\-^ , ■ ■ ■ , Q\„ such 
that Xi < ••• < A„, Ai G (— oo, Oi) U (fti, a2), and 
Xi G (ai_i,ai) U (aj,ai+i), for i = 2, . . . ,n. 

Set Uq = 0. If a generic point q is in the interior of the 
ellipsoid Q, then fxi > 0, so oq < /io < ai- In the same 



way, if a generic line i has a transverse intersection with 
the ellipsoid Q, then Ai > 0, so Ai G (ao,ai) U (01,02). 
The values (Iq = and Ai = are attained just when 
q (z Q and £ is tangent to Q, respectively. 

We denote by q = (/io, . . . ,/in) G M"+^, the Jacobi el- 
liptic coordinates of the point q = {xi, . . . , .x„+i). Carte- 
sian and elliptic coordinates are linked by relations 



Ut^j (aj -a-r)' 



(5) 



Elliptic coordinates define a coordinate system on each 
of the 2"+^ open orthants of the Euclidean space M"+^, 
but they become singular at the n -\- 1 coordinate hyper- 
planes, because the map g 1— )■ q is not one-to-one in any 
neighborhood of theses hyperplanes. 

A point is generic, in the sense of Theorem 2, if and 
only if it is outside all coordinate hyperplanes. From 
Eq. (5), we deduce that when the point q tends to the 
hyperplane Hj, some elliptic coordinate fXi tends to Oj. 
In fact, i = j OT i = j — 1, because of the inequalities 

Oi < Hi < flj+i. 

A line is generic, in the sense of Theorem 2, if and only 
if it is neither tangent to a singular confocal quadric nor 
contained in a nonsingular confocal quadric. 

If two lines obey the reflection law at a point q E Q, 
then both lines are tangent to the same confocal quadrics. 
This shows a tight relation between elliptic billiards and 
confocal quadrics: all lines of a billiard trajectory inside 
the ellipsoid Q are tangent to exactly n confocal quadrics 
Qai , • . • , Qa„ , which are called caustics of the trajectory. 
We will say that A — (Ai,...,A„) G M" is the caustic 
parameter of the trajectory. 

Definition 8 A billiard trajectory inside a nondegener- 
ate ellipsoid Q is nonsingular when it has n distinct non- 
singular caustics; that is, when its caustic parameter be- 
longs to the nonsingular caustic space 



A 



(Ai,...,A„) G 



< Ai < • • • < A„ 

Ai G (ai_i,Oi) U (aj,ai+i) 



We will only deal with nonsingular billiard trajectories 
along this paper. 

Remark 3 The set A has 2" connected components, be- 
ing each one associated to a different type of caustics. 
For instance, in the 2D case (n = 1) the two connected 
components correspond to ellipses and hyperbolas. 

Theorem 3 If a nonsingular billiard trajectory closes af- 
ter I bounces, all trajectories sharing the same caustics 
also close after I bounces. 

Poncelet and Darboux proved this theorem for ellipses 
and triaxial ellipsoids of M.^, respectively. Later on, it 
was generalized to any dimension in Ref. 14. 
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B. A dual transformation 

Let us present a map g : M ^i' M that takes place 
only for billiards inside ellipsoids. It exchanges the role 
of positions and velocities — so, it could be seen as a dual 
transformation — , and was introduced in Refs. 14 and 16. 
The following explicit formulae arc required to define g. 

Proposition 3 The billiard map f : M ^ M inside the 
ellipsoid (4) is expressed by ,f {q,p) = {q' ^p'), where 



{Dq,p') 
{Dp',p')' 
{Dq,p) 
{Dq,Dq)- 

Besides, the reversors r^f : M ^^ M are given by: 



q' = q + fi{q,p')p', fi{q,p')^-2 
p' =p + v{q,p)Dq, v{q,p) = -2- 



r{q,p) = {q,-p), q = q + ii{q,p)p, 
r{q,p) ^ {q,p), p= -p-v{q,p)Dq. 



Proof The formulae for the billiard map are well known. 
See, for instance, Ref. 16. Next, to obtain the formula 
for f we recall that q is the previous impact point, so 



q = q + tJ-{q,-p} (-p) = q + iJ'{q,p)p- 

Finally, the formula for f follows from p = —p'. 



We are ready to introduce the dual transformation g. 
We check that, in certain sense, g is the square root of /. 
In particular, g has the same symmetries and reversors 
as /. From our point of view, its most useful property is 
that it interchanges some symmetry sets. 

Proposition 4 Let g : M ~^ M , g {q,p) — {q,p), where: 



q = Cp' = Cp + iy{q,p)C-^q, 
p^-C-^q, 



with C = D ^'■^. Then, the following relations hold: 

a) .9^ = -/; and f og^ go f. 

b) g o So- ~ s^ o g for all a Cz S. 

c) f o g = —g of, g o f o g ~ f, and g o f o g ~ f. 

d) g (Fix (f^)) = Fix (f_^) = Fix (/ o f_„) for all cr G E. 
Proof 

a) We observe two relations: 



v{q,p) 



{Dq,p) 



{Dq,p') 



{Dq,Dq) {Dp',p') 

= -iJ-{q,p') = i^{q,p), 

,. _. ^ {Dq,p) ^ {Dq.p') 

^^'^^^^--^mp)-^TDqM 

= -v{q,p') = v{q,p). 



Thus vo g = jjLof and jio g — i/or, and then g^ = — /: 

ff' {q,P) - .9 {q,P) = {Cp + u {q,p) C-'q, -C-'q) 
= {-q-li{q,p')p',-p') = {-q',-p'). 

It is immediate that g is odd. Therefore, 



a ° (-.9^) = 9° f- 



f°9^ (-5^) ° .9 = -.' 

b) It is obvious. 

c) By using a) and definition of ; 



g{q,p)^{Cp + iy{q,p)C-'q,-C-'q) 

== {-q- i^{q,p)P,p) = ~f{q,p). 

As a consequence f o g = —g o f and besides: 

f = f o f = —g of = gofog^ 
f = fofof = fof = fo (—g ) 
= — f o g = g o r o g. 



r-| d) Because 5 is a diffeomorphism we have that: 



m e Fix (fa) ■<=^ g (m) = g {f„ (m)) = —fa {g (m)) 
<^=^ g (m) e Fix (f_^) . 

We have used that —fg- = f-^- Consequently, 

g (Fix (f,)) = Fix (f_,) = Fix (/ o f_,) . Q 

This proposition has a practical corollary. 

Corollary 1 The dual transformation g : M ^f M gives 
an explicit one-to-one correspondence between f„-SPOs 
and (/ o f_a)-SPOs of the same period. 

Proof Let m E Fix (f^) be a point of the phase space 
such that its orbit O, is Aperiodic. Thus, O is a Tct-SPO. 
Let us consider the orbit O by the point fh — g (ra) . 
First, 

f {m) = /' {g (m)) = g {f (m)) = g (m) = to. 

Thus, O is Z-periodic. Second, to — g (m) € Fix(f_cr), 
so O is a f_(j-SPO. Therefore, we have explicitly con- 
structed the correspondence g : O t^ O, which is one-to- 



one since g 



-/• 



D 



In subsequent sections we carry out some computations 
on SPOs only for reversors {fa- : ct G S}, since from this 
corollary, we deduce the analogous results for reversors 
{fa = f ofa -.a € E}. 
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C. Elliptic coordinates of SPTs 

The behaviour of elliptic coordinates along nonsingular 
billiard trajectories inside a nondegenerate ellipsoid, can 
be summarized as follows. If a particle obeys the billiard 
dynamics, describing a polygonal curve of R"+^ whose 
vertexes are on the ellipsoid, then its z-th elliptic coordi- 
nate fXi oscillates inside some interval li = [c2i,C2i+i] in 
such a way that its only critical points are attained when 
/ii G dli. In other words, the oscillation has amplitude 
C2i+i — C2i. This is a classical result that can be found, for 
instance, in Ref. 22. Let us state it as a formal theorem. 

In order to describe the intervals /,-, we set 



when they are expressed in elliptic coordinates. Besides, 
a billiard trajectory drawn in elliptic coordinates, has 
elastic reflections with the n-dimensional face {/ig = 0}, 
of the cuboid, but inner tangent contacts with the other 
2n + l faces. This behaviour can be observed in Table V, 
where several SPTs are drawn in Cartesian and elliptic 
coordinates. The cuboid is just a rectangle in those cases. 
Based on extensive numerical experiments, it has been 
conjectured^^ that 



2 < TO„ < • • • < mi < mo, 
but we are not aware of any proof. 



(6) 



{ci, . . . ,C2„+i} = {ai, . . . ,a„+i} U {Ai 



1 A„} . 



_i, and caustic pa- 



once ellipsoid parameters ai, . . . 
rameters Ai. . . . , A„, are fixed. 

If A e A, then ci, . . . ,C2n+i are pairwise distinct and 
positive, so we can assume that 

Co := < ci < ••• < C2„+i. 

Then li = [c2i,C2i+i], < i < n, arc the intervals that 
we were looking for. 

Theorem 4 Let q{t) be an arc-length parameterization 
of a nonsingular billiard trajectory inside the ellip- 
soid (4), sharing caustics Qai , • • • , Qa„ • Let q{t) — 
{lJ,o{t), . . . , ^j,n{t)) be the corresponding parameterization 
in elliptic coordinates. The following properties hold: 

fl) C2i < fJ^i (t) < C2i+i for all i G M. 

b) Functions fii (t) are smooth everywhere, except /ig (t) 
which is nonsmooth at impact points — that is, when 
q (i*) G Q — , in which case (Iq (t*+) = —Mo i^*^) 7^ 0. 

c) ///ii (t) is smooth att = t^, then 

H'i (t^) = -^=^ Hi (t^) G {C2i, C2j+i} . 

d) If the trajectory is periodic with length Lq, then q (t) 
is La-periodic and there exist some positive integers 
TTio, • . . ,m,„, called winding numbers, such that Hi (t) 
makes exactly mi complete oscillations (round trips) 
inside the interval li — [c2i,C2i+i] along one period 
< t < Lq. Besides, 

i. TTio is the period of the billiard trajectory. 

ii. mi is even when {c2i, C2i+i} H {ai, . . . , a„_|_i} 7^ 0. 

Hi. q {t) has period Lq/2 if and only if all winding 
numbers are even, but it never has period Lq/A. 

iv. Not all winding numbers can be multiples of four. 

Therefore, all billiard trajectories sharing the caustics 
Qai, • • • , Qa„ are contained in the (n + l)-dimcnsional 
cuboid 

Ca := /o X • • • X /„ = [0, ci] X • • • X [c2„, C2„+i] C M"+\ 



V. 2D CASE 

A. Caustics and elliptic coordinates 

We adapt the previous setting of billiards inside ellip- 
soids of IR"+-^ to the 2D case; that is, n = 1. Then the 
configuration space is an ellipse Q, which, in order to 
simplify the exposition, we write as 



q= {x,y) G 



X 

a 



y 



1 



a> b> 0. 



As we said in Subsec. IV A, any nonsingular billiard tra- 
jectory inside Q is tangent to one confocal caustic 



= 1 



\ = <q= (x,y) e 



a — A 6 — A 



where the caustic parameter A belongs to the nonsingular 
caustic space 



A^ EUH, 



E={0,b), H^(b,a). 



The names of the connected components of A come from 
the fact that Q\ is an confocal ellipse for A G -E and a 
confocal hyperbola for X d H. The singular cases A = 6 
and X = a correspond to the x-axis and y-axis, respec- 
tively. We say that the caustic type of a billiard trajectory 
is E or H, when its caustic is an ellipse or a hyperbola. 
We also distinguish between E-caustics and H-caustics. 

We denote by q = (e, h) the Jacobi elliptic coordinates 
of the point q — {x,y). That is, q belongs to the or- 
thogonal intersection of the confocal ellipse Qe and the 
confocal hyperbola Qh. We recall from Eq. (5) that 



x^ = 



(a — e) (a — h) 
a — b 



a — b 



Besides, 0<e<6</i<aif the point q is contained in 
the interior of Q, and e = at impact points q Cz Q. The 
y-axis is {h = a}, the x-axis is {e ~ b}U{h ~ b}, and the 
foci of the ellipse verify e — h = b, in elliptic coordinates. 
Likewise, the points (±x, ±y) have associated the same 
elliptic coordinates (e, ft.), so this system of coordinates 
do not distinguish among the four quadrants in M^. 
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FIG. 2. Geometric meaning of sides and vertexes of Cx 



We also know that all billiard trajectories sharing the 
caustic Q\ are contained in the rectangle 



C\ = [0, min (6, A)] x [max (6, A) , a] C 



(7) 



when they are expressed in terms of (e, h). Moreover, the 
coordinates e and h have a monotone behaviour except 
at the endpoints of the intervals that enclose them. The 
geometric meaning of the sides and vertexes of the rect- 
angle Cx is described in Fig. 2. In the figure, black, red 
and gray sides means impact with the ellipse Q, tangency 
with the caustic Qx, and crossing some coordinate axis, 
respectively. This color code is repeated in Table V. 



B. Rotation number and winding numbers 

We recall some concepts related to periodic trajectories 
of billiards inside ellipses. These results can be found, for 
instance, in Refs. 14 and 27. 

To begin with, we introduce the function p : A — > M 
given by the quotient of elliptic integrals 



P(A) 



/c 



niin(b, A) 




^(X-s)ib-s)ia-s) 



2r 

Jni 



max((),A) y(A_s)(b-.s)(a-s) 



It is called the rotation number and characterizes the 
caustic parameters that give rise to periodic trajectories. 
To be precise, the billiard trajectories sharing the caustic 
Qx are periodic if and only if 

p(A) = mi/2mo e Q 



TABLE I. Points (g,p) G Fix(r-) such that the line q + (p) is 
tangent to the caustic Qx- 

r q={x,y),p={u,v) 



R 
Rx 

Ry 

foRx 

foRy 

f ° Rxy 



rj.2 _ a{a.-X) 

g=(0,±^), 
g= (±Vo,0), u 
x^ = aX/b, y^ = 
2 ), „2 - 



a-b ' ^ 



:M(f,f)'M=\MA 



X 



a 



u — {a ~ A)/a, V — A/a 
' = A/6, v^ = {b-\)/b 
b-X,p= (±1,0) 
--b\/a,p= (0,±1) 



u 



-b ' 



V 



a-b' y 



x/ohjXp 



for some integers 2 < m\ < mo, which are the winding 
numbers. It turns out that mo is the period, mi/2 is the 
number of turns around the ellipse Qx for E-caustics, 
and mi is the number of crossings of the y-axis for H- 
caustics. Thus, mi is always even. Besides, all periodic 
trajectories with H-caustics have even period. 



C. Reversors and their symmetry sets 

We also change the notations for reversors. We de- 
note by R the reversor previously written as r. Be- 
sides, Rx, Ry, and R^y are the reversors obtained by 
composing R with the symmetries defined on the phase 
space and associated to the reflections {x,y) i->- {—x,y), 
{x,y) 1-^- {x,-y), and {x,y) <-)■ {-x,-y), respectively. 
Finally, f o R, foR^, foRy, and / o R^y denote the four 
reversors of the form {fa — fof„:a€'S}. We proved 
in Proposition 2 that the symmetry sets of Rxy and foR 
are empty sets, and following the same proposition we 
get: 



• Fix(i?) = {{q,p) €M:p€ N^Q}, 

• Fix(i?,) = {{q,p) eM:q= {0,±Vb)}, 
. Fix(i?^) ^ {{q,p) eM:q= (±^^,0)}, 

. Fix(/ o Rx) = {iq,p) GM:p^ (±1, 0)}, 

. Fix(/ o Ry) = {{q,p) e M : p = (0, ±1)}, and 

• Fix(/ o Rxy) = {iq,p) eM:Oeq+ {p)}. 

Here, M — {{q,p) E Q x §> : p directed outward Q at q} 
is the phase space. We write q = (x, y) and p — (u, v). 

D. Characterization of STs 

In order to find caustics Qx, whose tangent trajectories 
are periodic, we solve equation p (A) = mi/2mo, for some 
winding numbers 2 < mi < m-o. Hence, once we fix a 
caustic parameter A such that p (A) = mi/2mo £ Q, and 
a reversor r : M ^ M, it is natural to look for points 
{q,p) & Fix (r) such that q + (p) is tangent to Qx- If we 
can find those points, billiard trajectories associated to 
them are r-SPTs of period mo- 
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TABLE II. Points in STs that correspond to vertexes of C\ 



Vertex (e, ft) 



q* belongs to 



Rcversor 



Type 



(0,A) 


Qr\Qx 


R 


H 


(0,a) 


Q n {j/-axis} 


R. 


any 


(0,6) 


Q n { X-axis} 


Ry 


E 


(A, a) 


Qa n {y-axis} 


foR^ 


E 


(A,b) or (6, A) 


Q\ n {x-axis} 


foRy 


any 


(6, a) 


{(0,0)} 


foR.y 


H 



In the 2D case there are only 6 (= 2^ — 2) nonempty 
symmetry sets where one can look for STs; see previous 
subsection. We write the six formulae for (q, p) in Table I. 
The cases where p is directed inward Q at q are excluded 
from this table, because we require {q,p) € M. Let us 
check the first 3 cases in Table I. The formulae for /oi?^., 
/ o Ry and / o Rr^-y follow from their dual reversors Ry, 
R^ and R, respectively. See Corollary 1. Suffice it to 
perform the changes m^ 4-> x^ /a and v^ <-> y"^ /b in the 
first formulae. This has to do with the fact that the dual 
transformation g : M ^f M, g{q,p) = {q,p), p = (w, w) 
verifies that u — —xj^fa and v = —yjyh. 

R. We look for {q,p) e M such that p G N^Q and q+{p) 
is tangent to Q\. If p e NqQ, then there exists /i > 
such that p — ii{x/a,y/b). If, in addition, q + (p) is 
tangent to Q\, then Q\ is a hyperbola and q e QnQ\. 
Here, we have used that confocal ellipses and confocal 
hyperbolas are mutually orthogonal. Thus, the elliptic 
coordinates of g = {x,y) are (e, /i) = (0,A), and so 
x2 = a{a-X)/{a-b) and y^ ^ b{\ ^ b) / {a - b). 
Finally, we recall that p is a unit velocity: 



1 2,2 2 

1 = U + V = II 



r 

fe2 



ab 



^t■ = 



ab 
1' 



Rx- We look for a velocity p = (u, w) G § such that q+{p) 
has exactly one point on Q\ for q ~ (0, ±a/&). This 
condition turns into the uniqueness of solutions in t of 



where 



ar + 2/3t + 7 = 0, 



±Vbi 



A 6-A' 



&-A 



b-y 



or equivalently /32 = aj. This last equality together 
with u^ + v'^ = 1 gives u^ = A/a and u^ = (a — A) /a. 

Ry. It is obtained directly from the previous case by ex- 
changing the role of coordinates x and y. 

Next, we characterize nonsingular STs inside noncir- 
cular ellipses, as the trajectories passing, in elliptic coor- 
dinates, through some vertex of the rectangle C\. 

We only have to check that any vertex q^ e Ca, corre- 
sponds to some point q^, G q + (p), such that q + (p) is 



TABLE III. Feasible and forbidden reversors for each caustic 
type. 



Type 

E 
H 



Feasible 

Rx, Ry, f ° Rx, f ° Ry 
R, Rx, f ° Ry, f ° Rxy 



Forbidden 

R, f ° Rxy 
Ry, f ° Rx 



tangent to Qx and {q,p) G Fix(r) for some rcversor r. 
We list such correspondences in Table II. The point g^, 
is not unique, because elliptic coordinates do not distin- 
guish among the four quadrants in M^ . In its last column 
we describe the feasible reversors for each caustic type; 
see Subsec. VE for more detailed information. For in- 
stance, let us focus on the vertex q^, = (0,A). It is a 
vertex of C\ if and only if 5 < A < a, so the type of caus- 
tic is H. We note that qi, ^ Q since e = 0, and g^ G Q\ 
since h ~ X. Let p G § be the unique outward velocity 
such that <z* + (p) is tangent to Q\. Thenp G Nq^Q, since 
the confocal hyperbola Q\ has a orthogonal intersection 
at g* with the eUipse Q. Hence, {q*,p) G Fix(i?). The 
other cases are similar. We only observe that e 7^ in 
the last three cases, so g* ^ Q', but g^, is the middle point 
of two consecutive impact points. 



E. Forbidden reversors for each type of caustic 

Although there are 6 nonempty symmetry sets, once 
the caustic type (E or H) is fixed, only 4 of them give rise 
to SPTs of that type. The correspondences between caus- 
tic types and reversors whose symmetry sets are feasible 
or forbidden are listed in Table III. To check that all 
couples caustic type/feasible reversor take really place, 
we show some examples in Table V. 

We can prove Table III in three different ways. Firstly, 
by means of analytical arguments based on the formulae 
listed in Table I. Secondly, we could deduce it by using 
elliptic coordinates, as shown in Table II. Finally, we 
could write a geometric proof. The first way is the sim- 
plest one. If {q,p) G Fix(i?), then y^ == b{\ - b)/{a - b), 
so X > b and Qx is a hyperbola. If {q,p) G Fix(/ o R^y), 
then w^ ~ {X — b)/{a — b), so X > b and Qx is a hyperbola. 
If {<1,P) € Fix(i?y), then v'^ ^ {b - X)/b, so X < b and Qx 
is an ellipse. If {q,p) G Fix(/ o R^), then y^ — b — X, so 
X < b and Qx is an ellipse. This ends the proof. 



F. Characterization and classification of SPTs 

Next, we give a complete classification of nonsingular 
SPTs inside noncircular ellipses. We classify them by 
caustic type (E or H) and symmetry sets they intersect. 
We also characterize them as trajectories connecting, in 
elliptic coordinates, two vertexes of the rectangle Cx- 

Along this subsection, O denotes a r-SO through a 
point {q,p) G Fix(r), for some reversor r : M —)' M. We 
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TABLE IV. Symmetry properties of the billiard orbits 
fe:Pj) = F{<1,P) such that {q,p) £ Fix(r). 



r 


Symmetry properties 


R 


1-3 = 9j, P-j+i = -Pj 


R. 


<l-j = (-2;j,J/j), P-j = (itj+i,-uj+i) 


Ry 


Q'-j = fe>-yj)> p-j = (-'"j+i.t'j+i) 


foR^ 


9-0+1) = i-^j^y]), p-3 = ("j.-^'j) 


foRy 


Q-U+i) = (^j.-yj). P-J = (-'"j.^'j) 


foR^y 


'S'-O+i) = -<lj,P-3 =Pj 




set (gj,Pj) = fHq,p), qj = [xj.Vj), and p^ = iuj,Vj). 
The properties listed in Table IV are easily deduced using 
geometric arguments from the symmetry of the ellipse. 
For instance, let us explain the case of reversor Ry. If 
{q,p) e Fix(i?j,), then xq — ±y/a and j/o = 0. Let us 
assume that xq = ^/a. Then it is geometrically evident 
that q-j — {xj,—yj) and the billiard configuration is 
symmetric with respect to the x-axis; see Fig. 3. 

We note that all properties concerning velocities pj can 
be directly deduced from the properties of impact points 
qj, by using identity (2). 

Below, we suppose O is periodic as well, with winding 
numbers 2 < mi < tuq, mg being its period. We need a 
technical lemma on SPOs of even period. 

Lemma 2 Let us assume that the period is even: tuq — 
21, and the reversor is r — R^, so that xq — 0. 

a) // the caustic is an ellipse, qi ~ —Qq and pi — —po- 
Thus, qi-\-j — ~qj and pi^j = ^Pj , for all j G Z. 

b) If the caustic is a hyperbola, then qi — (—1) qq and 

(mq, —Vq) , if I odd and mi/2 even, 
pi = { {-uo,vo) , if I even, 

(— Uqj ^^o) 7 '/ ^ odd and mi/2 odd. 

Proof The hypothesis r — R^ in a), implies x^j = —Xj] 
see Table IV. Then, as toq — 21, we have xi = x^i = —xi, 
so xi = Q and qi e {go,-9o}- If qi = 90, then pi = po, 
since the line qi + (pi) must be tangent to the ellipse Q\ 
and the billiard trajectory turns around Q\ in a constant 
— clockwise or anticlockwise — direction. But the period 
cannot be smaller than m^ by hypothesis. Hence, qi — 
—go- Finally, the identity pi = ~po is necessary to keep 
the rotation direction around Q\. The formulae for a 
general integer j are obvious. 

We skip the proof for H-caustics in b): it is similar. D 

Following we prove that all SPOs inside a noncircu- 
lar ellipse are doubly SPOs and that there are 12 classes 
of SPOs, as claimed in the introduction. We recall that 
SPOs are classified by caustic type and intersected sym- 
metry sets. 



FIG. 3. A _R„-ST. 



Theorem 5 All SPOs inside a noncircular ellipse are 
doubly SPOs. SPTs inside noncircular ellipses are char- 
acterized as trajectories connecting, in elliptic coordi- 
nates, two different vertexes of their rectangle C\ . There 
are 12 classes of SPTs, listed in Table V. 

Proof Let O be a r-SPO through a point {q,p) G Fix (r) 
with winding numbers 2 < mi < ttiq, being ttiq its period. 
We must prove that there exists another reversor r ^ r 
and some index j £1, such that {qj,Pj) & Fix(r). 

First, let us consider the caustic type E. Then, accord- 
ing to Table III, we have only four possibilities: 



e {R^,foR,^,Ry,foRy}. 



(8) 



The easiest situation is when mo = 2fc + 1, because we 
know from the last item in Theorem 1 that 



{q,p)eFix{R, 



{qk+i,Pk+i) e Fix(/ o R^ 



and the same occurs with the couple of associated rever- 
sors Ry and f o Ry. This proves the first row of Table V. 
The formulae in the last column are deduced from the 
symmetry sets in Subsec. VC. 

Next, we assume that mo = 4fc. That is, we deal with 
the second row of Table V. Wc study each one of the 
four feasible reversors listed in Eq. (8). 



• li r = Rx, then go G {y-axis}, so xq 
mo = 4fc, we have y^k = y-k = yk ■ 



0. Thus, since 
-ysk. The last 

equality follows from Lemma 2a) for I = 2k and j — k. 

The above equalities imply y^^ = 2/fc = and thus we 

can take r — Ry and j G {fc, 3k}. 

• The case r = Ry is similar. Suffice it to exchange 
variables x and y. 

• The cases r E {f o R^, f o Ry} follow directly from the 
two previous ones by using Corollary 1. 

Finally, we assume that too = 4fc -I- 2, which corresponds 
to the third row of Table V. As before, the four feasible 
reversors listed in Eq. (8) are studied separately. 

• If r = Rx, then go G {y-axis}, so xq = 0. Thus, 
since too = Ak -\- 2, we have u^k+2 = ^-fc — u^+i — 
—U5k+2- The last equality follows from Lemma 2a) 
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ioT I ~ 2k + 1 and j = k + 1. These equalities imply 
U3fc+2 — Uk+i — and thus we can take r = f o Ry 
and j e {A: + l,3fc + 2}. 

• The case r = f o Ry follows directly from the case 
r = Rx by using Corollary 1. 

• The cases r € {Ry, f o Rx} are similar to the previous 
ones. Suffice it to exchange variables x and y. 

The proofs for the three rows about H-caustics follow the 
same lines, but using Lemma 2b). We skip the details. 

Once we know that all SPOs arc doubly SPOs, we 
deduce that, any doubly SPT connects two different ver- 
texes, namely, the ones corresponding to the couple of 
reversors of the doubly SPT. Indeed, we simply take into 
account the 1-to-l correspondence between reversors and 
vertexes of rectangles C\, as listed in Table II. Besides, 
there are 12 = 2 x 6 classes of SPTs, since there are two 
types of caustic (E and H) and any rectangle has four 
vertexes, and so, six couples of vertexes. 

Finally, we check that the 12 classes are realizable. 
That is, we ffnd one SPT of each class. This can be 
achieved by properly choosing the winding numbers. See 
Table V. ' D 



Therefore, there exist SPTs of even period mo = 21 > 
4, with only I + 1 distinct impact points on the ellipse. 
We have found similar SPTs inside ellipsoids of M.^ . See 
Fig. 1. 

Finally, we stress that any SPT has a dual version of 
the same period and caustic type, given by the dual trans- 
formation g; see Corollary 1. For instance, this duality is 
evident for the two SPTs with caustic type E that have 
common period ?no = 4fc. Suffice it to compare the two 
groups of associated formulae in the second row of Ta- 
ble V, and to recall that g exchanges the role of positions 
and velocities. In particular, the rectangular trajectory 
of that row is the dual version of the diamond-shaped 
one. The same applies to rows 1, 4 and 5. Instead, tra- 
jectories in rows 3 and 6 are dual of themselves. 



VI. 3D CASE 

The previous section sets the basis of this one. Roughly 
speaking, we apply similar ideas to find analogous results. 



A. Caustics and elliptic coordinates 



G. Minimal SPTs 

We end the study of SPTs inside ellipses by showing 
an SPT of each one of the 12 classes with the small- 
est period toq. We call minimal such SPTs. Therefore, 
mo E {3, 4, 6} for E-caustics, and mo G {4, 6} for H- 
caustics. Period mo = 2 is discarded, since two-periodic 
billiard trajectories are singular. The 12 minimal SPTs 
are drawn in Tabic V, both in Cartesian and elliptic co- 
ordinates. All of them connect two different vertexes of 
the rectangle C\, as claimed in Theorem 5. 

We realize that any minimal SPT of type E has a 
"twin" of type H when both are depicted in elliptic co- 
ordinates. This has to do with the fact that both tra- 
jectories have the same rotation number and connect the 
same vertexes. 

Any SPT, minimal or not, can be computed as follows. 
First, we fix a caustic type (E or H), a reversor r, a period 
mo (which must be even for H-caustics), and an even 
winding number mi, such that 2 < rrii < mg. We have 
listed in Table V all possible combinations. Then we 
find the caustic parameter through numerical inversion 
of the relation p (A) = mi/2mo. Of course, we solve that 
equation for A € i? or A G -ff, depending on the caustic 
type we are looking for. Finally, we take from Table I the 
starting point {q,p) G Fix (r) such that q+ (p) is tangent 
to Qx, and iterating the billiard map f : M ^ M, we 
get the whole SPO in the phase space M and the whole 
SPT in R2. 

Any i?-SPT is travelled twice in opposite directions, 
since it hits orthogonally the elhpse at points q G QCiQx. 
It can be observed in the last three rows of Table V. 



Let us recall some classical facts about caustics and 
elliptic coordinates. All the notations introduced in this 
subsection will be used later on in this section. 

We write the ellipsoid as 



(xi,X2,X3)GM-^^ + ^ + ^ = l^, (9) 

fli a2 as 



where < ai < a2 < a3. Any nonsingular billiard trajec- 
tory inside Q is tangent to two distinct confocal caustics 
Qai and Qx^ (with Ai < A2) of the family 



Qp = < (a;i,X2,a;3) G 



3 ^2 
p3 . Y^ -^j 



^-^ aj — u 



where /i ^ {oi, 02, 03}. If we set 

i;=(0,ai), iJi = (ai,a2), i?2 = (02,03), 

then Qp is an ellipsoid for fi Cz E, an one-sheet hyper- 
boloid if /i G -ffi, and a two-sheet hyperboloid when 
AiG i?2- 

Not all combinations of caustics exist. For instance, 
both caustics cannot be ellipsoids. The four possible 
combinations are denoted by EHl, HlHl, EII2, and 
H1H2. Hence, the caustic parameter A = (Ai,A2) be- 
longs to the nonsingular caustic space 

A={ExHi)U {Hi ® Hi) U{Ex H2) U [Hi x H2) , 

for Hi® Hi ^ {(Ai, A2) G -ffi X _ffi : Ai < A2}, in order 
to avoid the singular case Ai = A2. 
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TABLE V. Classification and symmetry properties of SPTs together to the minimal SPTs in Cartesian and elliptic coordinates. 
Black, red, grey, blue and green colors agree in both coordinates. The first three colors also agree with the ones displayed in 
Fig. 2. 



Type Period 



2k + 1 



E 



H 



4fc 



4fc + 2 



4/C + 2 
mi 



Ak 



4fc + 2 
mi 



odd 



Minimal SPTs in 
Cartesian coordinates 








Minimal SPTs in 
elliptic coordinates 








Couples of reversors 

and symmetry properties 

i?., foR, 

go = q2k+i G {y-axis} 

Pk+i II {a^-axis} 

Ry^ f ° Ry 

go = q2k+i G {a::-axis} 

Pfe+i II {j/-axis} 



Rx, Ry 

go = qik = — g2fc G {y-axis} 
gsfc = -gfc G {x-axis} 

f oR^, foRy 

Po = P4fc = -P2fc II {x-axis} 

Pik = -Pfe 1 1 {y-axis} 



Ry, / o i?^. 

go = g4fc+2 = — g2fc+i G {2;-axis} 
P3fc+2 = -Pfc+i II {x-axis} 

Rx, f ° Ry 

go = g4fc+2 = -g2fe+i G {y-axis} 
P3fc+2 = -Pfe+i II {j/-axis} 



Rx, f O Rxy 

go = g4fc+2 = -g2fc+i G {y-axis} 

gfc+i = — gfe, g3fc+2 ~ — gsfc+i 
xsk+i = Xk, yzk+i = —yk 

R, f O Ry 

go = ixo,yo) = g4fc+2 G Q n Qa 

g2*!+i = (a;o, — yo) € QnQ\ 

Pak+2 = -Pfe+i II {y-axis} 



^, fix 

go = g4fc = g2fe G {y-axis} 

Qk = {xk,yk) & Q f^Qx 

gsfc = {-Xk, yk) e QnQx 

foRy, foR^y 

go = g4fe = — gi, g2fe+i = — g2fc 
X2k = —xo,y2k = yo 
P3fc = —Pk 1 1 {y-axis} 

R, f ° fixy 
go = g4fe+2 = -g2fc+i G Q n Qa 

g3*!+2 = <lk ~ —Qk+l ~ —qsk+l 
Rx, f ° Ry 

go = g4fe+2 = — g2fe+i G {y-axis} 
P3fc+2 = -Pfc+1 II {y-axis} 



Let us write the three coordinates planes of M'^ as 

Ui = {(xi, X2, 2^3) eR^ -.xi^O}, le (1, 2, 3} . 

We know (see Theorem 2) that given any q ^ niUn2Un3, 
there exist an ellipsoid Qe, a one-sheet hyperboloid Qhi 
and a two-sheet hyperboloid Q/j^ such that 



g G Qe n Qh, n Qh^ 



(10) 



those being quadrics mutually orthogonal at q. There- 
fore, e < ai < hi < a2 < h2 < a^. Let q = (e, /ii, /12) be 



the elliptic coordinates oi q = (xi,X2,xs). The relation 
between Cartesian and elliptic coordinates (5) becomes 



(g,; - e) (fli - hi) [aj - /12) 
(oi - flj) (a, - flfc) 



(11) 



for any {i, j, k\ — {\, 2, 3}. To understand what happens 
at nongeneric points, we define the three singular caustics 
as 



n. 



le{i,2,3}. 
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We also recall that the focal ellipse and focal hyperbola 
shared by the family of confocal quadrics {Q^ : /i G M} 
are 



= 1 CHi, 



£ = <^ (0,a;2,X3) G 



H = Uxi,0,X3)e 



a2 — 0,1 03 — fli 



U C Ho. 



as — 0,2 0-2 — ai 

These notations are useful because of the following 
lemma, which is a known extension of the first item in 
Theorem 2. 

Lemma 3 If q ^ £ WH, there exists q = (e, hi, /12) such 
that e < ai < hi < 02 < h2 < 03, e 7^ hi, hi 7^ /i2; and 
property (10) holds, including its orthogonal character. 

Thus, we can express in elliptic coordinates all points 
along nonsingular billiard trajectories. The trajectory is 
contained inside Q, so e > 0. Impact points q € Q cor- 
respond to e = 0. Tangency points q G Q\ . correspond 
to e = Xj E E, hi ^ Xj £ Hi or /i2 — Xj E H2 de- 
pending on the type of the caustic Q\^ . Crossing points 
q E Hi have the following elliptic coordinates: Hi is 
{e = ai} U {hi = ai}, 112 is {hi = a2} U {/12 = 02}, and 
lis is {/l2 = 03}- 

Indeed, one could agree that £ = {e = ai = hi} and 
Ji = {hi — 02 = h2}, although we do not need it, because 
all nonsingular billiard trajectories sharing the caustics 
Qai and (5^2 are contained, in elliptic coordinates, inside 
the cuboid 



Cx = 



[0,Ai] X [ai, A2] X [02, a3] , for EHl, 

io,ai] X [Ai,A2] X [02,03], for mm, 

io,Ai] X [01,02] X [A2,03] , for EH2, 

[0,ai] X [Ai,02] X [A2,03] , for mH2. 



(12) 



We note that e 7^ ft-i 7^ /12 for all q = (e, /ii,/i2) G Cx, 
since A = (Ai, A2) G A. Hence, all points belonging to Cx 
verify Eqs. (10) and (11). 



B. Frequency map and winding numbers 

The following results related to periodic trajectories of 
billiards inside ellipsoids can be found in Rcf. 27. 

There exists a map uj : A — > M^, called frequency map 
and defined by means of six hyperelliptic integrals, that 
characterizes all caustic parameters A = (Ai, A2) that give 
rise to periodic trajectories. (Its explicit formula is given 
in Ref. 27.) To be precise, the billiard trajectories sharing 
the caustics Qxi and Qx2 are periodic if and only if 



u (Ai, A2) = (mi, m2) /2mo e 



(13) 



for some positive integers toq, mi, and m2. These are 
the winding numbers introduced in item d) of Theorem 4. 
They describe how trajectories fold in M"^. First, mo is 
the period. Second, mi is the number of Hi-crossings and 
m.2 is twice the number of turns around the xi-axis for 



TABLE VI. The symmetry sets of the 16 reversors of the 
billiard inside a triaxial ellipsoid Q C K"*. 



Fix (r) 



Fix(/ o r) 



R 


peNgQ 





Ri 


q £ Si and p G NqSi 


Urn — Un = 0,U^ = 1 


-tt'TTin 


Xjn = Xn ^ U, Xi = ai 


{q+ {p))±,{xi-&xis} 


R123 





0€q + {p) 



EHl-caustics; mi is twice the number of turns around the 
X3-axis and m2 is the number of H3-crossings for EH2- 
caustics; mi is the number of tangential touches with 
each one-sheet hyperboloid caustic and m,2 is twice the 
number of turns around the xi-axis for HlHl-caustics; 
whereas m-i is the number of H2-crossings and 7712 is the 
number of H3-crossings for HlH2-caustics. Only EHl- 
caustics and EH2-caustics have trajectories of odd pe- 
riod. 



C. Reversors and their symmetry sets 

Let CT;, cTmn, 17123 ■ ^^ ^ 11^ '^ be the cCmXn-specular re- 
flection, the xi-axial reflection, and the central reflection, 
respectively. Here, {l,m,n} ~ {1,2,3}. Two samples 
are CT3 (xi,a;2,a;3) = (a;i,a;2, -2:3) and cri3 (a;i,X2,a;3) = 
(-xi, 2:2, -2:3). Let si,s.mn,si23 : M -> M be the 
symmetries given by si{q,p) = {aiq,aip), Smn{q,p) = 
{crmnq,crmnP), and Si23{q,p) = {cri23q,ai23P)- We de- 
note by R the rcversor written as f in Sec. HI. Likewise, 

Rl ^ SiO R, Rjnn = cr„i„ o i?, and i?i23 = S123 o -R- 

Remark 4 Henceforth, we shall adopt the following in- 
dex notation for the sake of brevity. If I, m, and n ap- 
pear without an explicit explanation of their ranges, it 
must be understood that they represent any choice such 
that {I, m, n{ — {1, 2, 3} . This convention is used, for in- 
stance, in Tables VI and VII. Moreover, in the second to 
fifth rows of Table VII, index I has a distinguished mean- 
ing from indexes m and n; and the formulae contained in 
those rows are invariant under permutation of m and n. 
This has to do with the fact that Umn = fnm- 

The three planar sections of the ellipsoid are 

Si^QfMVu ?G{1,2,3}. 

The symmetry sets of the 16 — 2'^ reversors are listed 
in Table VI, which follows directly from Proposition 2. 
We note that only the symmetry sets of i?i23 and f o R 
are empty. We usually write coordinates q = (xi,X2,x^) 
in the configuration space Q, and p = [ui, U2, M3) in the 
velocity space §^. 



D. Characterization of STs 

We want to characterize nonsingular STs inside non- 
degenerate ellipsoids, as trajectories passing, in elliptic 
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coordinates, through some vertex of the cuboid Cx. In 
order to accomplish this goal, we must prove both impli- 
cations. In Proposition 5 we check that any vertex of C\ 
is visited by some (indeed, eight) STs. In Proposition 6 
we see that any ST passes through some vertex of C\. 

Proposition 5 Let A = (Ai,A2) G A &e fixed. Given 
any vertex q^ = (e, /ii,/i2) of Cx, let q^, € Qe(^ Qhi n 
Qh2 ■ Then there exists some reversor r, and some point 
{q,p) G Fix (r) such that g* G l+ip) and the line q+{p) is 
tangent to Qai o,nd Q\2- If e = 0, then q = qi,; otherwise 
q^, is the middle point of the consecutive impact points q 
and q. 

We clarify some aspects of this result before its proof: 

• We follow the index notation explained in Remark 4. 

• We know that two symmetry sets only intersect at 
points {q,p) € M such that the hne q+{p) is contained 
in some coordinate plane; see Proposition 2. But we 
are dealing with nonsingular trajectories. Therefore, 
we can associate just one reversor to each vertex. 

• Correspondence q^ h^ (QtP) is l-to-8; sec Table VII. 
This has to do with two facts. First, elliptic coordi- 
nates do not distinguish among the eight octants in 
M^. Moreover, unit velocities p directed inward Q at 
impact points q are excluded from Table VII; other- 
wise, {q,p) ^ M. The bilhard trajectory associated to 
any of those eight points {q,p) € M is r-symmetric and 
passes, in elliptic coordinates, through vertex q^,. 

Proof To begin with, we know from Lemma 3 that Qe^ 
Qhi and Qh2 are mutually orthogonal at q*. There arc 
several kinds of vertexes. The first distinction is e = or 
e ^ 0; that is, qi, & Q or q^, ^ Q. The second distinction 
is the number of coordinate planes through q^,. It turns 
out that there are six kinds of vertexes. We study each 
kind separately. 

a) If the point g* belongs to Q but is outside all three co- 
ordinate planes, then q* = (0, Ai, A2). From Eq. (12), 
the type is H1H2. Thus, Qai is a one-sheet hyper- 
boloid, Q\2 is a two-sheet hyperboloid, and Q, Qai 
and Q\2 are pairwise orthogonal at q ~ q^,. Thus, if 
p is the outward unit normal velocity to Q at q, then 
{q,p) € Fix (i?) and the line q + (p) is tangent to Qx^ 
and Qx2- Finally, let us check the formulae for q and 
p given in Table VII. 

The formula for q follows directly from Eq. (11). Thus, 
there are exactly eight choices for g. li p G NqQ and 
p points outward Q at q, then there exists /i > such 
that ui = ^xi/ai. And using that p G §^ wc get that 



1=M^ 



/i = 



aia2a3 



A1A2 



In summary, only _R-STs with 111112-caustics can take 
place for this vertex. 



b) If the point qi, belongs to Q and is contained in only 
one coordinate plane, say 11;, then q^, = (0, /ii,/i2) 
with {/ii,/i2} = {a;,Aj} for some j G {1,2}. In that 
case, q^: e Q n IIj n Qx^ = SiCi Qx^ ■ We do not get any 
restriction on the type of the caustics, but only some 
restrictions on the indexes; sec Eq. (12). Namely, 

• / e {2, 3} and j = 2 for type EHl; 

• le {2,3} for type HlHl; 

• / e {1, 2} and j = 2 for type EH2; and 

• l=j = 2 or {l,j) ^ (3, 1) for type H1H2. 

We observe that any {q,p) E M such that q ^ qi, and 
p G NqSi is contained in Fix(_R;). Thus, the formula 
for q follows again from relations (11), and we find four 
choices for q. Next, wc look for an outward velocity 
p — (ui, U2, U3) G §^ such that pUTqSi and q + {p) is 
tangent to both Qx and Qx^. with {.?', fc} = {1,2}. 

Condition pUTqSi reads as x„Mm/a„ = XmUn/am- In 
other words, we look for some i^ G M such that 

Then, tangency with Qx is immediate, since Qx is 
orthogonal to the planar section Si — Q O Hi at q. It 
only remains to impose the tangency with Qx^. ■ Thus, 
the quadratic equation 



{xi+tUl) {X„i+tu„i) (Xn+tUn) 



ai — \k 



^m. '^h 



a„ — \h 



1 



must have a unique solution in i, which is equivalent 
to the vanishing of its discriminant: /?^ — 0:7, where 



1 



a; — Afe a„j — Afe a„ — Aj. aj — A^ 
1 



O-l ^m 2 2 , ^l ^n 2 2 \ 2 



'^m'^r, 






a^nOniX- I3)v 



^k {a,n — Afc)(a„ — Afe) ' 



- 1 



O-m ^ Afe Cln — Xk 

_ Xk{Xj — Xk) 

(Om — Xk){an — Xk) 

Here, we have used that uf — 1 — uf^ — u'^. After 
some calculations we find v^ — Xk/Xjabc, and from 
here we get the desired formulae in Table VII. We 
must impose in turn, that p is an outward velocity, so 

< {Dq,p) = I' (unxlJa^ + a-mX^/an) , 

and v > 0. Thus, once we fix q, there are two choices 
for ui, and exactly eight choices for {q,p). 

We finish this item by stressing the consequences of 
the restrictions on the index I. They mean that the re- 
versor i?i cannot take place for types EHl, HlHl, and 
H1H2, whereas R^ is forbidden for type EH2. This in- 
formation is displayed in Table IX. 
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TABLE VII. Relations among vertexes q* of the cuboid C\, types of caustics, reversors r, and points q-,, £ q + (p) such that 
q+{p) is tangent to Qx^ and Q\2, with {q,p) G Fix (r). Notation for indexes {l,m,n} = {1, 2, 3} and {j,k} = {1, 2} is described 
in the text. 



Vertex q* = (e, /ii,/i2) 


Type 


q* belongs to 


r 


q = (xi.Xg.Xs), P= (wi,M2,M3) 




(0,Ai,A2) 


H1H2 


QnQx^nQx^ 


R 


2 _ ai (ai — Ai) (ai — A2) _ xi 
{ai — a„i) {ai — a„) ' a; \ 


1 A1A2 



= 0, {hi,h2} = {ai,\j} all Q n n, n Qa, Ri 



e^O, {hi,h2} ^ {arn,a„} jjj^jjj^ Q H 11^ H n„ R„ 



xi = 0, X 



2 

5 '''m 



^m \^m ^j ) 



O'n yO,n Aj J 



0.1 — Afc 



A, 



ai V aia2a[i\j 

2 n 2 -^1-^2 . . X 
Xi = a;, a;^ = a;„ = 0, u, = , sign (ui) 



■ sign (x; 



(a„ - Ai) (a„ - A2) 2 _ (fi" ~ -^i) (ffln - ^2) 



(^m yf^rn (^n ) 



Ojn [ Ctn Q,rj 



Wj = 1^ Um = Un = 0, :r^ = , Sign (ui) = sign (xi) 

2 _ {am — Ai) (am — A2) 2 _ (<^n — Ai) (ttn — A2) 

"^ m, 1 •*-' n. 



{e,/ii,/i2} = {a;, Ai, A2} 



not 
HlHl 



H; n Qai n Qa2 / o ft 



{e,/ii,/i2} = {a„,a„, Aj} 



all n^ n n„ n Qa, / o i?„ 



2 "'"T, ^*J 

-uz = 0, U^ = , U 

^l ^ 0,1 — Afc, [XmfXn) = 



aman-^k 



aiXj 



\'-^m 1 "'n j 



(01,02,03) 



H1H2 ninn2nn3 f^R^r, u} = ]"' ' ^'\^''' ' ^\ 

(ai — am) [ai — a„) 



xi = 



010203 



A1A2 



-ui 



If Qi, belongs to Q and is contained in two coordinate 
planes, say 11^ and n„, then q^, — (0, ft.i,ft.2) with 
{^1,^2} = {drmO-n}- This prevents HlHl-caustics; 
see Eq. (12). Furthermore, we get again some restric- 
tions on the indexes. Namely, 

• 1 e {m, n} for type EHl; 

• 3 e {m, n} for type EH2; and 

• 1 ^ {m, n} for type H1H2. 

Besides, q* is one of the two vertexes of the ellipsoid 
on the x/-axis. Hence, ((j*,p) € Fix(_Rmn) for any 
{q*,p) € M. Finally, let us check that the formulae 
given for p — (ui, 1*2,^3) in Table VII are correct. 

We look for an outward velocity p G S^ such that the 
line q + (p) is tangent to Q\^, for k = 1,2. Thus, the 
quadratic equation 



{xi+tUif {x,n+tu,nf (x^ + tu,iY 



ai 



At 



— At 



a„ — A/( 



1 



must have a unique solution in t, which is equivalent 
to the vanishing of its discriminant: /3^ = aj, where 



P = 



ai - Afe 


a-m 


— Afe a„ — Afe 




1 


( 

1 + 
V 


a-a-ai 2 , an - 


-am 


a-n - Afc 


ai - Xk am- 


-Xk 


XlUl 








ai- >^k' 








t2 


1 — 


Afe 





ai 



ai 



After some calculations we obtain that 

On (flm - Afe) wf + {am - ««) AfeU^j == Afe («„ - Afe) . 

for k = 1,2. These two equations along with condition 
uf + u^ + M^j = 1 form a system of three linear equa- 
tions in the unknowns uf , u^, and u^. A simple com- 
putation shows that the formulae given in Table VII 
are the unique solution of this system. 

Finally, we impose that p is an outward velocity, so 

< {Dq,p) == xiui/ai, 

and xiui > 0. Thus, once fixed the vertex q, there are 
two choices for Um and two others for u„. Therefore, 
there are exactly eight choices for {q,p). 

We finish again by stressing the consequences of the 
restrictions on the indexes m and n. They mean that 
the reversor R23 cannot take place for type EHl, R12 
is forbidden for type EH2, both R13 and R12 are not 
allowed for type H1H2, and all reversors of the form 
Rmn cannot be found for type HlHl. This informa- 
tion is also displayed in Table IX. 

d) If qi, does not belong to Q but it is in all coordinate 
planes, then g^ = (0,0,0) = and q^ = (01,02,03). 
In that case the type is H1H2; see Eq. (12). We look 
for lines g* -I- (p) = (p) that are tangent to the one- 
sheet hyperboloid Qai and to the two-sheet hyper- 
boloid Qx^- Since g^, is the center of both quadrics, 
p = (mi,M2,U3) must be an asymptotic direction of 
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them, so 






ul , 


ul 


«l 


ai ~ Afe a2 


-A, ' 


as - -^fc 


These equations along 


: with u^ 



fce{i,2} 



4 + u„ = 1 form a 
system of three linear equations in the unknowns uf, 
u^, and u^. A simple computation shows that the 
formulae given in Table VII are the unique solution of 
this system. Thus, there are eight choices for p. 

The line q^, + (p) = + (p) intersects the ellipsoid Q 
at some points q — fxp and q ~ — /ip where /i > 0. So, 
xi = iiui. And using that g G Q we get that 



1 = M^ 



/i = 



aia2a3 
A1A2 



Finally, we realize that {q,p) G Fix(i?i23) and q^, is 
the middle point of q and q. 

In summary, only (/ o _Rx23)-STs with IIlH2-caustics 
can take place for this vertex. 

If (7vf does not belong to Q but it is in two coordinate 
planes, say ![„ and n„, then q^ = (e, /ii,/i2) with 



{e,hi,h2} = {a„ 



,Xj}- This case is the "dual" of 



case b), so it is associated to the reversors / o Rmn- 
The computations do not involve new ideas. We skip 
them. We just stress that there are no restriction on 
the type of the caustics, but only some restrictions on 
the indexes; see Eq. (12). Namely, 

• 1 e {m, n} for types EHl, HlHl, and H1H2; and 

• 3 e {m, n} for types EH2. 

They mean that the reversor /oi?23 cannot take place 
for types EHl, HlHl, and II1II2, whereas / o R12 is 
forbidden for type EII2. 

f) If q^, docs not belong to Q and is contained in only 
one coordinate plane, say 11/, then q^, = {e,hi,h2) 
with {/ii,/i2} = {aZjAi,A2}. In that case, q^, Cz Tli D 
Q\i n Q\2- This case is the "dual" of case c), so it is 
associated to the reversors / o i?;. We also skip the 
computations. 

We only remark that no HlHl-caustics can be asso- 
ciated to this case. Besides, we get again some re- 
strictions on the indexes; see Eq. (12). Namely, / 7^ 1 
for type EHl; Z 7^ 3 for type EH2; and / = 1 for type 
H1H2. Consequently, /oi?i cannot take place for type 
EHl, /oi?3 is forbidden for type EH2, both /oi?2 and 
foR^ are not allowed for type H1H2, and all reversors 
of the form f o Ri cannot be found for type HlHl. D 

Remark 5 Opposed vertexes of a given cuboid Cx pro- 
vide points on dual symmetry sets. At the same time, the 
formulae for {q,p) of reversors foRi, foRmn and foRi23 
follow from the formulae of their dual reversors Rmn : Ri 
and R, respectively, by the change uf O xf/ui. This does 



TABLE VIII. Symmetry properties of the arc-length param- 
eterization q{t) of r-STs such that g(0) = g« and </'(0) = p. 

r Symmetry properties 



R 
Ri 

foRi 

f oRm 



q{~t) = q{t) 
q{-t) =(Jiq{t) 
q{—t) = o-mnq{t) 
q{-t) =(Jiq{t) 



H-t) 



.q{t) 



foRi23 q{-t) = -q{t) 



not come as a surprise; see Corollary 1. Sufhce it to re- 
alize that the dual transformation g : M — > M, g {q,p) ~ 
(q,p), verifies that p = —C^^q, C^ = diag (ai, a2, a3). 

All implications in the proof of Proposition 5 can be 
reversed. Hence, we can move along Table VII in both 
directions: from left to right, and from right to left. 

Nevertheless, we prefer to prove Proposition 6 follow- 
ing a reasoning whose generalization to an arbitrary di- 
mension is straightforward. 

Proposition 6 Let O be a nonsingular r-SO through a 
point {q,p) G Fix (r) for some reversor r : M ^- M. Let 
A = (Ai,A2) & A be its caustic parameter. Then there 
exists a vertex q* = (e, /ii,/i2) of the cuboid C\ and a 
point q^ G Qe n Qhi H Qh2 such that q^, <E q+ {p) . 

Proof First, we consider the reversors of the form i?, Ri 
or Rran- Let g : M ^ M? , q{t) = {xi (i) ,X2 (i) ,^3 (i)), 
be any of the arc-length parameterizations of the bil- 
liard trajectory through the impact point qi, := q with 
unit velocity p, determined by the conditions q{0) — q^, 
and g'(0) = p. Clearly, it is smooth except at impact 
points. Besides, it has the symmetry properties listed in 
Table VIII, which can be deduced from comments after 
Remark 2 and elementary geometric arguments. 

Let q : M ^ Cx, q (i) = (e (t) , hi (t) , /12 (t)), be the cor- 
responding parameterization in elliptic coordinates. We 
know that components e(t), hi (t), and ft.2 (t), oscillate 
in some intervals in such a way that their only critical 
points are attained at the extremes of these intervals. 
The cuboid C\ is the product of these three intervals. 
Besides, the functions hi (t) and /i2 (t) are smooth every- 
where. 

The above symmetries of g (i) imply that q (— i) — q{t), 
since elliptic coordinates do not distinguish among the 
eight octants in R^. Therefore, 

h'l (0) = 0, h'2 (0) = 0. 

In addition, e (0) = 0, since we have taken qi, = q <E Q. 
Hence, q* — (e, /ii,/i2) := q (0) = (0, ft-i (0) , /ii (0)) is a 
vertex of C\] see item c) of Theorem 4. 

It remains to consider reversors of the form f o Ri / o 
Rmn and / o i?i23. Let g : R — >• R^ be the arc-length 
parameterization of the billiard trajectory that begins at 
the middle point q^, := (q + q) /2 with unit velocity p. 
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TABLE IX. Forbidden reversors for each type of caustics. 



Type 



Forbidden reversors 



EHl 
EH2 

HlHl 

H1H2 



R, fo Ri23, Ri, f o R23, R2S, foRi 
R, fo R123, R3, f o R12, R12, foR3 
R, Ri, R12, Ria, R23, f o ^123 
f o Ri, f o R2, f o R3, f o R23 
Ri, f o R23, Rl3, f o R2, R12, / o i?3 



We recall that q is the previous impact point, that is, 
Qniq+{p))^{q,q}. 

The symmetry properties of these parameterizations 
are also listed in Table VIII. For instance, if r = / o 
i?i23, then q^, — (0,0,0), which implies q{—t) ~ —q{t). 
Then, we apply exactly the same argument as before, 
with just one difference. The function e (i) is smooth at 
t — 0, because q{Q) — qi, ^ Q and e (0) ^ 0. Hence, we 
also get that e' (0) — 0, so q* — (e, /ii,/i2) := q (0) = 
(e (0) , hi (0) , hi (0)) is a vertex of Cx. 

Finally, the intersection property q* e Qe n Qhi H Qh2 
follows from Lemma 3. D 

Corollary 2 Nonsingular STs inside triaxial ellipsoids 
of M^ are characterized as trajectories passing, in elliptic 
coordinates, through some vertex of their cuboid. 



E. Forbidden reversors for each type of caustics 

Next, we emphasize that, although there are 14 non- 
empty symmetry sets, once fixed the caustic type, some 
of them cannot take place. To be more precise, there 
are ten forbidden symmetry sets for HlHl-caustics, but 
only six otherwise. These results have been obtained 
along the proof of Proposition 5, and we organize them 
in Table IX. Forbidden reversors appear in dual couples 
— couples whose symmetry sets are interchanged by the 
dual transformation g. 

A glance at cuboid (12) and Table VII shows that, 
two different vertexes of the same cuboid cannot be as- 
sociated to the same reversor, but for HlHl-caustics. 
This agrees with the previous discussion, where we found 
8 — 14—6 feasible symmetry sets for caustics of type EHl, 
EH2, and H1H2, but only 4 = 14- 10 for HlHl-caustics. 



F. Characterization and classification of SPTs 

Let us characterize and classify SPTs inside triaxial 
ellipsoids of M?. We classify them by the caustic type 
and the couple of vertexes of the cuboid they connect. 

Theorem 6 Nonsingular SPTs inside triaxial ellipsoids 
of M.^ are characterized as trajectories connecting, in el- 
liptic coordinates, two different vertexes of their cuboid. 
All nonsingular SPOs are doubly SPOs, but a few with 



HlHl-caustics. Besides, there are exactly 112 classes of 
nonsingular SPTs, listed in Tables X-XIII. 

Proof Let O be a nonsingular r-SPO through a point 

{q,p) G Fix(r), for some reversor r : M ^ M, with 
winding numbers TOo,mi,TO2. Let A — (Ai,A2) G A be 
its caustic parameter. Let Lq be its length. 

Let q : M — >■ M'^ be the arc-length parameterization of 
the billiard trajectory that begins at the distinguished 
point q^, introduced in Proposition 6. 

Let q : M — > Ca be the corresponding parameterization 
in elliptic coordinates. Then q (t) is even — see the proof 
of Proposition 6 — and Lo-periodic — see Theorem 4. 

The key trick is as simple as to realize that 

q {Lo/2 -t) = q{t- Lo/2) = q{t + Lq/2) , Vi G M. 

Hence, q (t) is even with respect to Lq/2, and so 

q, = (e*, hi, hi) := q (0) , q. = {e' ,h\, h^) := q (Lo/2) 

are vertexes of the cuboid C\. The proof for q^, was al- 
ready explained in Proposition 6; the proof for q, is equal. 
Now, taking into account the interpretation of the 
winding numbers as, the number of complete oscillations 
of each elliptic coordinate when the arc-length parameter 
t moves from to Lq — see item d) of Theorem 4 — , we 
consider two cases: 

• If some winding number is odd, then q, 7^ q*, because 

e* = e* ^ mo G 2Z, h* = h'^ ^ rui e 2Z. 

• If all winding numbers are even, then q (t) has period 
Lo/2 and q, — q^, — see item d) of Theorem 4 — , in 
which case we repeat the same strategy to find a second 
vertex qo = {e°,h1, /12) '■— ^ (-^o/4) such that 

e* = e° 4^ mo G 4Z, h* = h° -^ m^ <E 4Z. 

Then, qo 7^ q* — see item d) of Theorem 4. 

Therefore, the trajectory connects two different vertexes 
in both cases. Next, we prove that any trajectory q{t) 
connecting, in elliptic coordinates, two vertexes is an 
SPT. Thus, we assume that q^ = q (i*) and q, = q (t.) 
are two different vertexes with t* < t,. From Propo- 
sition 5 we know that this trajectory is r^-symmetric 
and r,-symmetric for some reversors r* and r,. The 
case r* = r, is not excluded. Repeating the reason- 
ing in Proposition 6, we deduce that q (t) is symmet- 
ric with respect to t = i^, and t ~ t,. In particular, 
q (2t* + t) = q i-t) and q {2t, + i) = q (^0 for all t G K. 
Set T = 2{t,-t^). Then 

q{t + T) = q (2t, + t - 2t^) = q (2t* - t) = q (t) 

for all t G M. Hence, q (t) is T-pcriodic, and so, q [t] is 
periodic with period T or 2T; see item d) of Theorem 4. 
This proves the characterization of SPTs. 
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TABLE X. Classification of SPTs for caustic type EHl. No- 
tation for kinds "e", "o", "f", and "t" is described in the text. 
For each kind of winding numbers we hst its minimal repre- 
sentative, and its four couples of reversors, whose order inside 
the couple is irrelevant. 



TABLE XL Analogous of Table X for caustic type EH2. 



(mo,mi,m2) 
kind/minimal 



Couples of reversors 



(t,t,t) 
(10,6,2) 



{Ra,foRi2},{R2,.f°Ri3} 
{Ri3,foR2},{Ri2j°Rs} 



(t,t,f) 
(10,6,4) 



{R3joRr,},{R2,foRr2} 
{Ri3,foR3},{Ri2joR2} 



(t,f,t) 
(6,4,2) 



{R3,foR2},{R2,foR3} 

{Ri3,foRi2}, {Ri2,foRr3} 



(o,e,c) 
(5,4,2) 



{R3joRs},{R2,foR2} 

{Rl3,f'=-Rl3}ARl2,foRi2} 



(f,t,t) 

(8,6,2) 



{R3,Rl2},{f°R3,.f°Rl3} 
{R2,Rl3},{f°R2j°Rl3} 



(f,t,f) 

(8,6,4) 
(f,f,t) 
(8,4,2) 



{R3,Rl3},{f°R3,f°Rl3} 
{R2,Rl2},{f°R2,f°Rl2} 
{R3,R2},{f°RiJ°R2} 
{Rl3,Rl2}, {/oi?i3,/oi?i2} 



Once we have established that any SPT connects two 
different vertexes of its cuboid, it is easy to deduce that, 
all SPOs inside triaxial ellipsoids of R^ whose type of 
caustics are EHl, EH2 or H1H2, must be doubly SPOs. 
Suffice it to recall that each vertex of the cuboid is asso- 
ciated to a different reversor for those three types. 

The number 112 comes from 4 x 28, since there are 4 
types of caustics and any 3-dimensional cuboid has eight 
vertexes, and so 28 = (8 x 7) /2 couples of vertexes. 

We must check that none of the 112 classes of SPTs 
is fictitious. To do it, we present an algorithm that pro- 
vides a minimal SPT of each class. Minimal means that 
it has the smallest possible period. The key step of 
the algorithm is to properly choose the winding num- 
bers (TOo,mi,TO2). In accordance with the previous dis- 
cussion, we distinguish four kinds of winding number. 
Namely, even: "e" , odd: "o" , multiple of four: "f " , and 
even but not multiple of four: "t" . For instance, all wind- 
ing numbers of kind (t,t,t) connect opposite vertexes of 
their cuboids. The kind (f, f, f) never takes place — see 
item d) of Theorem 4. 

Let us explain the algorithm by using an example. We 
want to connect opposite vertexes q^, = (e*,/i^,ft.2) — 
(0, A2,a3) and q, ~ (e*,/!*,/!*) = (Ai, 01,02) of the 
cuboid 

Cx = [0,Ai] X [ai,A2] X [02,03], 

which corresponds to caustic type EHl. 

This problem about vertexes is equivalent to find an 
SPT of type EHl that is simultaneously _R3-symmetric 
and (/ o i?i2)-symmetric, because the reversors associ- 
ated to q* and q, are R3 and / o R12, respectively (cL 
Table VH) . As we said above, winding numbers must be 
of kind (t,t,t). Then, taking for granted conjecture (6), 
the minimal choice is (?71o,toi, 7712) = (10,6,2). Finally, 



{mo, mi, 7712) 
kind/minimal 



Couples of reversors 



(t,t,t) 
(10,6,2) 



{Ri,.foR23},{R2,foRr3} 

{Rl3,foR2},{R23joRl} 



(t,t,f) 

(10,6,4) 



{Ri,foR2},{R2,foR,} 

{Rl3,.f°R23}, {R23,foRr3} 



(t,f,t) 

(6,4,2) 



{Ri,.foRr3},{R2,foR23} 

{Rl3joRl},{R23joR2} 



(o,e,e) 
(5,4,2) 
(f,t,t) 
(8,6,2) 



{Ri,.foRi},{R2,foR2} 
{R13J0R13}, {R23J0R23} 

{Rl,R23},{f°Rl,f°R23} 
{R2,Rl3},{f°R2j°Rl3} 



(f,t,f) 

(8,6,4) 



{Ri,R2},{.foRi,foR2} 

{Rl3,R23}, {foRi3,foR23} 



(f,f,t) 

(8,4,2) 



{Ri,Ri3},{foRi,foRr3} 

{R2,R23},{f°R2j°R23} 



TABLE XII. Analogous of Table X for caustic type H1H2. 



(mo, mi, 7712) 
kind/minimal 



Couples of reversors 



(t,t,t) 
(10,6,2) 


{R,foRi23},{R2,f°Rl3} 

{R3j°Rl2}, {R23J0R1} 


(t,t,f) 
(10,6,4) 


{R,foRr2},{R2,foR3} 

{R3,foRi23}, {R23,foRl3} 


(t,f,t) 
(6,4,2) 


{R,foRis},{R3,foRi} 
{R2J0R123}, {R23J0R12} 


(t,f,f) 
(10,8,4) 


{R,foRi},{R2,foRi2} 

{R3,foRr3}, {R23,foRl23} 


(f,t,t) 
(8,6,2) 


{R,R23},{f°Rl,.f°Rl23} 
{i?2,i?3}, {f°Rl2joRl3} 


(f,t,f) 
(8,6,4) 


{R,R2}, {f Rr3, f R123} 
{R3,R23}, {foRi,foRi2} 


(f,f,t) 
(8,4,2) 


{R,R3},{.foRi2,foRi23} 
{R2,R23}, {/oi?i,/oi?i3} 



if we solve Eq. (13) for A = (Ai,A2) G E x Hi, with 
(toq, TOi, 7712) = (10,6,2), and draw the billiard trajec- 
tory through the point {q,p), given by the formulae of 
the second row of Table VH for I — 3, then we obtain 
a trajectory of period 10, type EHl, i?3-symmetric, and 
(/ o i?i2)-symmetric. 

All winding numbers of kind (t,t,t) connect oppo- 
site vertexes of their cuboids. Hence, their SPTs of 
type EHl can display the following double symme- 
tries: {i?3,/o^i2}, {-R2,/o^i3}, {■Ri3,/o-R2}, or 
{^12, / ° ^3}- They are obtained by looking at Table VH 
the relations between vertexes and reversors. 

Other kinds of winding numbers and other caustic 
types can be studied in a completely analogous way. And 
so, we complete Tables X-XHL D 

The caustic type HlHl has the following peculiarity: 
both caustics Qai and Qx^ are 1-sheet hyperboloids. We 
will denote Qx^ as the outer hyperboloid and Qx2 as the 
inner one, since Ai < A2. Note that two vertexes of the 
form (e, Ai, ^,2) and (e, A2, /12) are associated to the same 
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TABLE XIII. Analogous of Table X for caustic type HlHl. 
See comments after the proof of Theorem 6 for the meaning 

of (I). 



(mo, mi, 7712) 

kind/minimal 

(t,t,t) 

(10.6,2) 



Couples of reversors 

(i?3|/oi?l2), {R2\foRl3) 
{foRl2\R3), ifoRl3\R2) 



(t,t,f) 

(10,6,4) 



(i?2|/oi?l2), (i?3|/°i?13) 

if o Ri2\R2), if ° RnlRs) 



(t,f,t) 

(6,4,2) 



iR2,foRr,\),i\R2,foRi3) 
iR3joRi2\),i\R:i,f°Ri2) 



(t,f,f) 

(10,8,4) 



(i?2,/oi?i2|), (|i?2,/oi?i2) 
(i?3,/oi?i3|), (|i?3,/°i?13) 



(f,t,t) 

(8,6,2) 
(e,o,e) 
(4,3,2) 



iR2\R3), if o Rvilf o R^2) 

iRiR2), ifoRi2foRi3) 

(-R2I-R2), if o Ri2\f o R12) 
iRslRs), if o Rvilf o Rrs) 



(f,f,t) 
(8,4,2) 



iR2,R3\),if°Rl2,foRl3\) 
i\R2,R3),i\f°Rl2,foR,s) 



reversor; see Table VII. In particular, SPTs connecting 
vertexes of this form arc not doubly symmetric. Likewise, 
different connections can give rise to the same couple of 
reversors. We still classify those connections as different 
classes, because they have different geometries. Symbols 
(i?2,i?3|), (i?3|i?2), iR2\R3), and (|i?2,i?3) denote four 
classes with caustic type HlHl and the same couple of 
reversors: {i?2,-R3}- They are depicted in Table XVI. 
But the points g*, associated to each reversor (Table VII) 
are on the outer (respectively, inner) hyperboloid when 
the reversor is written to the left (respectively, right) of 
symbol |. This notation is used in Table XIII. 

We have seen that, once fixed the caustic type, there 
exists a tight relation between the symmetry sets associ- 
ated to an SPT and the kind of winding numbers. A sim- 
ilar result for standard-like maps was obtained by Kook 
and Meiss.^^ Compare Tables X~XIII with their Table I, 
where they list the 36 classes of SPOs for the 4D sym- 
plectic Froeschle map. Our 4D symplectic billiard map 
has more classes because a new factor enters in the clas- 
sification: the caustic type. 



G. Gallery of minimal SPTs 

We have 112 classes of SPTs listed in Tables X-XIII, 
so we tackle the task of finding a minimal representative 
in each class. They have periods toq <E {4, 5, 6, 8, 10}. 

The algorithm for the caustic type EHl is: 1) Choose 
one of the seven minimal winding numbers (toq, mi, TO2) 
in Table X; 2) Find caustic parameters Xi G E and A2 G 
Hi such that Eq. (13) holds; 3) Choose one of the four 
couples of reversors {r, r} in the corresponding row of 
Table X; 4) Get a point iq,p) from Table VII using the 
reversor r (or r); and 5) Draw the doubly SPT through 
q with velocity p. 



Only step 2) is problematic, because it requires the 
inversion of the frequency map. The main obstacle is 
that Eq. (13) might not have solution in E x Hi for some 
of the minimal winding numbers at hand. Nevertheless, 
it is known^^ that oj\exHi is a diffeomorphism and 



hm i^iE X Hi) = Uu}i,u}2) 

ai— >0+ 



< ^2 < CJi < i } . 



Consequently, Eq. (13) has a unique solution in E x Hi 
provided that ellipsoid (9) is fiat enough. 

The other three caustic types can be dealt with the 
same algorithm. There is just one difference among them. 
Namely, we must consider different shapes of ellipsoids to 
ensure that Eq. (13) has a solution. More precisely, we 
must consider almost flat ellipsoids (ai small) for caustic 
types EHl and HlHl, and almost "segments" (both oi 
and a2 small) for caustic types EH2 and H1H2. These 
results can be found in Ref. 27. This explains why some 
ellipsoids in Tables XIV~XVII are so flat. We have chosen 
not too extreme ellipsoids whenever it has been possible. 

We have found a minimal SPT of each class. Some of 
them are displayed in Tables XIV-XVII. Arrows in their 
column headings point in the direction of increasing xi, 
X2 or x^. Each line of Tables XIV, XV and XVII shows 
four different perspectives of an SPT. The 3D image cor- 
responds to an isometric view of the first octant. The 
three projected planes are viewed from the positive miss- 
ing axis. We only display 3D views in Table XVI. Bil- 
liard trajectories are depicted in red. Green and yellow 
lines represent intersections of Q with Hl-caustics and 
H2-caustics, respectively. The E-caustic is also shown in 
cases EHl and EH2. 

Period four takes only place for four classes of SPTs of 
type HlHl. Period five takes only place for eight classes 
of SPTs, half of type EHl and half of type EH2. We 
display half of these minimal SPTs in Table XIV. Their 
announced symmetries can be verified by observing their 
projections. We can check that trajectories of caustic 
type HlHl give one turn around the xi-axis and have 
three tangential touches with the outer (respectively, in- 
ner) 1-sheet hyperboloid, trajectories of type EHl give 
one turn around the xi-axis and cross four times the 
plane Hi, whereas trajectories of type EH2 give two turns 
around the Xa-axis and cross twice the plane H2. This 
is consistent with the geometric interpretation given in 
Subsec. VI B of winding numbers. 

We show in Table XV the minimal representatives of 
the classes listed in the third row of Table XII. Since 
all of them have (6, 4, 2) as winding numbers, they are 
6-pcriodic, cross four times the plane H2, and cross twice 
the plane H^. 

We draw several SPTs of type HlHl in Table XVI. 
They show the difference among the six classes (-R2 1-^2)7 

(i?3|i?3), (i?2,i?3|), iR2\R3), iR3\R2), and (|i?2,i?3)- We 

can locate those classes in the last rows of Table XIII. 
One can observe, for instance, that the SPT of class 
(i?2, -R3I) has two impacts on the intersection Qr\Il2r\Q\^ 
and two more on Q H3 H Qai • On the contrary, the 
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TABLE XIV. Two minimal SPTs of period 4 and four minimal SPTs of period 5. The "Data" column includes consecutively: 
caustic type, winding numbers (mo, mi, 7712), ellipsoid parameters a^ and ai — we have set 03 = 1 — , caustic parameters A2 
and Ai, and all reversors whose symmetry sets intersect the corresponding SPOs. 



3D {xi :t,X2 ■■\,X3 -.y/) 



Plane Hi {x2 :t, x^, 



Plane 112 {x\ :t, 2:3 :^) 



Plane Il-j, (ii :t, X2 



Data 





i 



HlHl 

(4,3,2) 

0.8 

0.13 

0.648376 

0.130077 

i?2 





r^,. 















HlHl 

(4,3,2) 

0.8 

0.13 

0.648376 

0.130077 

EH2 

(5,4,2) 

0.3969 

0.2 

0.762965 

0.199523 

Ri 

foRi 

EH2 

(5,4,2) 

0.3969 

0.2 

0.762965 

0.199523 

R23 
/0-R23 

EHl 
(5,4,2) 

0.49 

0.25 
0.260266 
0.231635 

R12 

foRi2 

EHl 
(5,4,2) 

0.49 

0.25 
0.260266 
0.231635 

R3 
foR, 
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TABLE XV. The four classes of SPTs with minimal winding numbers (mo, mi, mi) — (6,4,2) and HlH2-caustics. "Data" as 
in Table XIV. 



3D {xi ■.t,X2 ■■\,X3 :>/) 



Plane Hi {x2 '-^jXa :— >) Plane 112 {xi :t, a^s :<— ) 



Plane lis {xi :t, X2 



Data 





f#r 




H1H2 

(6,4,2) 

0.45 

0.13 

0.967756 

0.133273 

R 
foRia 






H1H2 

(6,4,2) 

0.45 

0.13 
0.967756 
0.133273 

i?2 
/ O i?i23 






H1H2 

(6,4,2) 

0.45 

0.13 
0.967756 
0.133273 

Ra 







H1H2 

(6,4,2) 

0.45 

0.13 
0.967756 
0.133273 

^23 
f O Rl2 



SPT of class (|i?2,i?3) has two impacts on the intersec- 
tion Q n 112 n Qx2 ^^^d two more on Q n Us n Qx^ . Since 
Ai < A2, Qxi and Qx^ are the outer and inner one-sheet 
hyperboloids, respectively. 

Finally, in Table XVII we present half of the 6-periodic 
minimal SPTs with an ellipsoidal caustic. They corre- 
spond to the third row of Tables X-XI. 

Some comments about these trajectories are in order: 

• The 4-periodic trajectories in Table XIV are the sim- 
plest examples of nonplanar periodic trajectories. Be- 
sides, they are one of the few simply SPTs. Both tra- 
jectories have two points on the same symmetry set. 
However, each point is associated to a different caus- 
tic: inner or outer. That property is easier to see in 
the first one. Hence, they are of class (i?2|^2) and 
(/ ° Ri^lf ° ^13) in the notation used in Table XIII. 

• The 5-periodic trajectories in Table XIV are the sim- 
plest examples of nonplanar periodic trajectories with 
odd period. They have one point on a symmetry set 



and another on the associated symmetry set. We recall 
that only SPTs with odd period have this property; see 
item c) of Theorem 1. 



• Many projections onto the horizontal plane Hi 
like 2D SPTs (cf. Table V). 



look 



• Any i?-SPT is travelled twice in opposite directions, 
since it hits orthogonally the ellipsoid at some point 
q e Q n Qxi n Qx2 ■ Therefore, there exist SPTs of 
even period toq — 21 > 6 with only I + 1 distinct im- 
pact points on the ellipsoid, although all of them are 
of caustic type H1H2. We show a 6-periodic sample 
(the simplest one) in the first row of Table XV. 



VII. THE GENERAL CASE 

We have characterized STs and SPTs in terms of the 
vertexes of some rectangles (for the 2D case) and cuboids 
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TABLE XVI. Six SPTs (3D view) for HlHl-caustics. They show different classes of SPTs corresponding to reversors R2 and/ or 
_R3. "Data" as in Table XIV. 

3D {xi ■■t,X2 ■\,X3 -.y/) I Data | 3D {x\ :t, a^2 :\, 3:^3 '-i/) \ Data | 3D (xi :t, 3^2 :\,a:3 '-^Z) \ Data 





HlHl 

(4,3,2) 

0.8 

0.13 

0.648376 

0.130077 

{R2\R2) 



HlHl 

(4,3,2) 

0.8 

0.13 

0.648376 

0.130077 

[R-AR'i) 





HlHl 
(8,4,2) 

0.95 

0.05 
0.056134 
0.457414 
{R2,Rs\) 



HlHl 
(8,4,2) 

0.95 

0.05 
0.056134 
0.457414 

(|i?2,i?3) 





HlHl 

(8,6,2) 

0.95 

0.05 
0.050041 
0.229595 
{R2\Rs) 



HlHl 

(8,6,2) 

0.95 

0.05 
0.050041 
0.229595 
{Ra\R2) 



TABLE XVII. Four minimal SPTs with winding numbers (mo, mi, 7712) = (6,4,2) and an ellipsoidal caustic. "Data" as in 
Table XIV. 

3D (xi ■■\,X2 '■\,X3 :■/') Plane Hi {x2 :t, 2:3 :— >) Plane H2 (xi :t, a^a :<— ) Plane Hs {xi :t, 3^2 :— >) Orbit data 





A 



EH2 

(6,4,2) 

0.45 

0.13 
0.962896 
0.126968 

R23 
foR2 






EH2 

(6,4,2) 

0.45 

0.13 

0.962896 

0.126968 

R2 
f ° R23 





EHl 

(6,4,2) 

0.8 

0.13 

0.403278 

0.126231 

R12 
f o Ria 





EHl 

(6,4,2) 

0.8 

0.13 

0.403278 

0.126231 

R2 

foR-i 
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(for the 3D case). We have proved that there are ex- 
actly 12 and 112 classes of SPTs in these low-dimensional 
cases. See Theorems 5 and 6. 

Next, we generalize both characterizations and the fi- 
nal classification to nondcgenerate ellipsoids of M"+^. 

We say that two SPTs are of the same class when they 
have the same type of caustics and they connect the same 
couple of vertexes of their respective cuboids. All SPTs 
inside the same class are associated to the same rever- 
sors. For instance, looking at Table II, we see that all 
SPTs whose type of caustic is E and that connect the 
two left (respectively, right) vertexes of the rectangle (7) 
are associated to the reversors R^ and Ry (respectively, 
f o R^ and /o Ry). 

We have the following result. 

Theorem 7 Nonsingular STs (respectively, nonsingular 
SPTs) inside nondegenetate ellipsoids ofM"'^^ are char- 
acterized as trajectories passing, in elliptic coordinates, 
through some vertex (respectively, two different vertexes) 
of their cuboid. All nonsingular SPOs are doubly SPOs, 
but a few with caustics of repeated types. There are 
2^" ^2"+^ — l) classes of nonsingular SPTs. 

This theorem is obtained by means of small refine- 
ments of the techniques used in this article, although 
some checks become rather cumbersome. 

The number 2^" (2"+-'^ — l) has a simple explanation. 
There are 2" types of caustics — see Remark 3 — , and 
any [n + l)-dimensional cuboid has 2"+^ vertexes, and 
so 2" (2"+^ — 1) couples of vertexes. 

Once fixed the ellipsoid and the caustic type, all 
2" (2"+^ — l) couples of vertexes are realizable. That is, 
each couple is connected by some SPT. Of course, one 
should choose suitable winding numbers for each couple. 
That choice is guided by the following observation. 

Let C\ be a cuboid such that its billiard trajecto- 
ries are periodic with winding numbers mo, . . . , ?Ti„. Let 
q — (/xq, . . . , An) and q = {jlo, . . . , /i„) be any couple of 
different vertexes of Cx connected by an SPT. If some 
winding number is odd, then fli = fii <^ mi Cz 2Z. If 
all winding numbers are even, then /ti = /l^ <=> jrii G 4Z. 
Here, i is any integer index such that < i < n. 



VIM. CONCLUSIONS AND FUTURE WORK 

We have shown that the billiard map associated to a 
convex symmetric hypersurface Q C K""*"^ is reversible. 
Even more, it admits 2"+^ factorizations as a compo- 
sition of two involutions. Therefore, its SPOs can be 
classified by the symmetry sets they intersect. We have 
carried out this classification for nondegenerate ellipsoids 
of M"+^. The characterization of STs and SPTs in terms 
of their elliptic coordinates has been the key tool. 

The existence of SPOs is useful for several reasons. We 
indicate just two. 

First, it reduces the numerical computations needed 
to find POs. Concretely, we can halve the dimension of 



the problem, by restricting the search of periodic points 
on the n-dimensional symmetry sets. (This issue will 
be greatly appreciated when dealing with perturbed el- 
lipsoids, for which no frequency map exists, since the 
integrability is lost.) Second, an SPO of a completely 
integrable reversible map, persists under symmetric per- 
turbations, when the symmetry sets are transverse to the 
Liouville invariant tori of the map at some point of the 
SPO. 

These two facts are crucial to study billiards inside 
symmetrically perturbed ellipsoids, which is our next 
goal. To be more precise, we plan to generalize the results 
about the break-up of resonant invariant curves for bil- 
liards inside perturbed circles^^ and perturbed ellipses^^. 
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